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Abstract 

This  thesis  investigates  MX  survivability  when  a  layered 
or  terminal  defense  system  is  deployed  with  various  numbers 
of  multiple  protective  shelters  (MPS) .  The  layered  defense 
system  defends  the  MX  with  an  exoatmospheric  layer  which  is 
augmented  by  an  endoatmo spheric  layer.  The  exoatmospheric 
layer  protects  the  MX  with  longwave  infrared  (LWIR)  guided 
interceptors  which  must  directly  impact  an  incoming  RV  at 
approximately  300,000  feet  altitude  to  destroy  it.  The  endo- 
atmospheric  layer  consists  of  a  terminal  BMD  system  known  as 
Low  Altitude  Defense  (LoAD)  which  defends  the  MX  with  three 
hypersonic,  nuclear  armed  interceptor  missiles.  The  terminal 
defense  system  consists  of  either  one  or  two  LoAD  systems. 

This  research  effort  determines  the  most  cost  effective  defense 
system,  and  draws  conclusions  on  these  systems  based  upon 
quantitative  and  qualitative  (ex.,  political)  considerations. 
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MX  SURVIVABILITY:  PASSIVE  AND  ACTIVE  DEFENSE 


I.  Background 


General  Topic 

According  to  a  statement  by  former  Secretary  of  Defense 
Brown  in  October  1980,  America's  land-based  ballistic  missile 
force  may  now  be  vulnerable  to  Soviet  ICBM  attack  (Ref  23:16). 
Two  causes  cited  are  the  rapidly  increasing  numbers  of  Soviet 
reentry  vehicles  (RVs)  and  technological  advances  allowing 
RV  accuracy  to  increase  to  unpredicted  levels  in  this  decade. 

In  1976,  the  U.S.  initiated  full  scale  development  of  a 
more  secure  basing  mode  consisting  of  MX  missiles  in  multiple 
protective  shelters  (MPS).  The  strategy  is  to  base  200  mis¬ 
siles  in  4600  horizontal  shelters,  with  one  missile  located 
in  each  group  of  23  shelters  (see  Figure  1).  After  hardening 
to  a  certain  level,  this  obscuration  of  a  missile's  precise 
location  appears  to  be  the  only  viable  passive  technique  left 
to  increase  U.S.  ICBM  survivability  (Ref  20:15).  The  probabil¬ 
ity  of  a  single  enemy  RV  destroying  one  MX  missile  cannot 
exceed  one  in  23  as  long  as  the  location  of  the  missile  within 
the  cluster  of  23  shelters  remains  unknown  (i.e.,  preservation 
of  location  uncertainty  is  maintained).  Although  President 
Reagan  has  proposed  that  the  MX  be  initially  based  in  existing 
silos,  this  proposal  has  yet  to  be  ratified  by  the  Congress. 
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Figure  1.  MX  MPS  Basing  Plan  (Ref  15) 
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Furthermore,  the  final  MX  basing  plan  has  not  been  determined, 
and  MPS  is  certainly  a  viable  option  which  needs  to  be  studied. 

As  advances  in  Soviet  technology  improve  the  effective¬ 
ness  of  each  RV,  the  probability  of  one  RV  successfully  destroy¬ 
ing  one  MX  shelter  (an  easily  locatable  target)  will  approach 
one.  At  this  point,  4600  RVs  could  destroy  all  4600  shelters 
and  the  200  MX  missiles.  Due  to  this  threat,  the  Army  (mili¬ 
tary  department  responsible  for  ICBM  defense)  is  currently 
studying  various  types  of  ballistic  missile  defense  (BMD) 
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systems.  This  study  will  investigate  MX  survivability  as 
a  function  of  several  possible  defense  systems. 

General  Situation 

Since  the  strategic  capability  of  the  two  major  military 
powers  has  drawn  equal  in  the  past  decade,  the  U.S.  strategy 
for  nuclear  deterrence  has  been  one  of  approximate  equality 
of  strategic  offensive  capability  (i.e.,  essential  equiva¬ 
lence)  (Ref  1:2).  Defensive  forces  have  not  been  a  part  of 
U.S.  Strategic  Policy  since  the  signing  of  the  Anti-Ballistic 
Missile  (ABM)  Treaty  with  the  Soviet  Union  in  1972  and  the 
abandonment  of  the  U.S.  operational  Safeguard  ABM  site  in 
North  Dakota.  The  U.S.  considered  defensive  weapon  systems 
that  might  threaten  Soviet  deterrent  capabilities  as  destabi¬ 
lizing  (Ref  1:2)  and  adopted  the  idea  in  the  1970s  that 
modernization  of  offensive  forces  and  BMD  development/deploy¬ 
ment  were  mutually  exclusive  policies  (Ref  20:12).  It  was 
felt  that  U.S.  strategic  deterrent  goals  could  be  met  with 
offensive  weapons,  exclusively. 

In  recognition  of  the  growing  strategic  imbalance  and 
vulnerability  of  our  ICBMs,  a  change  in  the  U.S.  nuclear 
targeting  policy  was  formalized  in  1980  (Ref  8:65-66).  It 
is  now  a  "countervailing"  strategy  that  allows  attack  of 
military  targets  (counterforce  targets),  industrial  targets/ 
population  centers  (countervalue  targets),  or  a  combination 
of  these  two. 
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In  order  to  retain  the  countervailing  strategy  (with 
its  implied  assured  destruction  and  damage  limiting  roles) 
and  still  promote  arms  reduction,  the  U.S.  strategic  force 
structure  must  be  robust  enough  to  withstand  technological 
surprise,  be  somewhat  insensitive  to  arms -control  "cheating," 
respond  economically  to  threat  growth,  and  finally,  promote 
crisis  stability.  It  does  not  seem  possible  to  achieve  these 
goals  with  strictly  an  offensive  force  structure  (Ref  8:3  and 
20:15). 

It  has  been  suggested  by  several  reliable  sources  and 
verified  by  initial  studies  that  a  properly  configured  stra¬ 
tegic  force  of  offensive  missiles  and  a  ballistic  missile 
defense  would  permit  continued  deterrence,  allow  meeting  the 
goals  stated  above,  and  still  have  significant  economic 
advantages  over  an  all-offensive  force.  Additionally,  this 
force  could  achieve  the  above  regardless  of  Soviet  behavior/ 
response  (Ref  8:3  and  20:11). 

For  example,  if  the  U.S.  deployed  two  ABMs  per  MX  cluster 
(i.e.,  one  MX  missile  and  23  shelters),  Soviet  strategy  must 
target  at  least  three  RVs  at  each  shelter  (two  to  defeat  the 
BMD  and  one  to  destroy  the  MX)  to  insure  destruction  of  the 
one  missile  in  an  MX  cluster  --  a  total  of  69  RVs.  Therefore, 
Soviet  destruction  of  the  planned  U.S.  MX  configuration  of 
200  MX  clusters  would  require  a  total  of  13,800  RVs.  With 
this  type  of  leverage,  it  is  possible  to  envision  a  reduction 
in  the  total  number  of  shelters  needed  for  a  preset  level  of 
ICBM  survivability. 


Due  to  the  impact  of  MX  on  ecology,  the  economy,  and 
the  defense  posture,  additional  study  of  MX  survivability  with 
and  without  various  defense  elements  is  beneficial  if  for  no 
other  reason  than  to  clearly  elucidate  all  possible  alterna¬ 
tives  within  technological  reach.  An  extensive  search  of 
unclassified  literature  (for  instance,  Refs  14,  17,  27)  gives 
no  indication  that  the  layered  BMD  problem  described  has  been 
researched  in  total  using  simulation  (see  Ref  22  for  a  simu¬ 
lation  of  the  low  altitude  component  only).  All  discovered 
references  use  mathematical  models  and  very  narrowly  defined 
Ballistic  Missile  Defense  systems.  The  problem,  as  defined 
in  this  paper,  lends  itself  very  well  to  the  use  of  simulation 
because  of  the  relatively  undefined  nature  of  this  postulated 
defense  system.  A  general  analysis  of  an  undefined  BMD  system, 
using  simulation,  will  be  a  helpful  first  step  toward  investi¬ 
gating  the  capabilities  of  a  particular  system  if  one  is 
defined  in  more  detail  later. 

Problem  Statement 

This  thesis  addresses  MX  survivability  with  and  without 
employment  of  a  layered  BMD  system,  as  well  as  the  cost 
effectiveness  of  survivability  alternatives.  A  layered  BMD 
system  consists  of  two  elements,  (1)  an  exoatmospheric 
component,  and  (2)  an  endoatmospheric  component.  The  endo- 
atmospheric  component  will  be  the  Low  Altitude  Defense  (LoAD) 
system  as  defined  by  James  Moore  (Ref  22). 


The  effect  of  the  layered  defense  on  MX  survivability  in 
many  varied  situations  will  be  examined  and  analyzed,  as  will 
the  effect  of  different  numbers  of  MX  shelters  per  complex. 
This  study  should  quantify  the  cost  effectiveness  and  surviva 
bility  of  defense  alternatives  as  the  MX  is  defended  by 
different  configurations  against  a  postulated  Soviet  attack. 
The  overall  results  should  delineate  the  effectiveness  and 
relative  costs  of  (1)  various  active  defenses  and  (2)  various 
numbers  of  MX  shelters  when  confronted  by  a  Soviet  threat. 


Overall  Objective 

The  objective  of  this  study  is  to  develop  a  methodology, 
using  simulation,  which  quantifies  MX  survivability  for 
various  combinations  of  passive  and  active  defense. 


Specific  Goals 

Four  specific  goals  have  been  established  for  this  study 

1.  Construct  a  computer  model  which  will  compute 
MX  survivability  under  certain  MPS  and  defense 
configurations,  and  a  predetermined  Soviet  threat. 

2.  Explore  the  tradeoff  relationship  between  the 
number  of  MX  shelters  per  complex  and  a  layered 
defense  system. 

3.  Explore  the  tradeoff  relationship  between  the 
number  of  MX  shelters  per  complex  and  terminal 
defense  systems  alone. 

4.  Analyze  the  results  and  comment  on  the  cost 
effectiveness  of  possible  solutions  to  the  U.S. 

ICBM  vulnerability  dilemma,  while  taking  into 
account  qualitative  factors  such  as  political 
implications . 
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Scope 

When  constructing  the  models,  it  will  be  essential  to 
calculate  the  probability  of  kill  (PK)  of  an  attacking  Soviet 
RV  and  the  PK  of  the  U.S.  interceptors  when  launched  at  an 
incoming  RV.  To  accomplish  this,  four  data  sets  (RV,  MPS, 
terminal  defense,  and  layered  defense)  which  contain  the 
information  required  to  perform  the  necessary  calculations 
must  be  developed.  Once  values  are  selected  for  the  required 
parameters,  the  effectiveness  of  various  combinations  of 
passive  (MPS)  and  active  (terminal  or  layered)  defense  tech¬ 
niques  will  be  examined. 

The  research  will  begin  with  a  description  of  the  Soviet 
threat  and  the  U.S.  defense  techniques  utilized  in  this 
analysis.  A  description  of  the  Soviet  threat  will  include 
such  factors  as  the  number  of  incoming  RVs,  their  yield,  and 
CEP;  while  a  description  of  the  U.S.  defenses  will  include 
such  systems  as  the  MPS,  the  terminal  defense,  the  exoatmos- 
pheric  defense,  and  the  layered  defense. 

The  following  alternatives  will  be  examined: 

1.  Given  one  terminal  BMD  system  per  MX  complex  and 
various  numbers  of  MPS  per  MX  complex,  determine 
the  number  of  exoatmospheric  interceptors  per 
complex  necessary  to  obtain  a  predetermined  level 
of  survivability. 

2.  Given  one  terminal  BMD  system  per  complex  and  no 
exoatmospheric  defense  layer,  determine  the  num¬ 
ber  of  MPS  per  complex  necessary  to  obtain  a 
predetermined  level  of  survivability. 
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3.  Given  two  terminal  BMD  systems  per  complex  and 
no  exoatmospheric  defense  layer,  determine  the 
number  of  MPS  per  complex  necessary  to  obtain  a 
predetermined  level  of  survivability. 

During  the  analysis  phase,  the  most  significant  variables 
of  each  data  set  can  be  identified  while  the  effectiveness 
of  the  passive  and  active  defense  systems  is  determined. 

Limitations 

In  order  to  facilitate  this  study,  several  assumptions 
have  been  made.  Some  of  these  are  based  on  actual  facts, 
while  others  are  made  to  assist  in  the  modeling  process. 

The  number  of  interceptors  per  terminal  defense  unit  is 
assumed  to  be  three,  to  be  consistent  with  the  recent  study 
of  James  Moore  (Ref  22).  It  is  also  assumed  that  the  U.S. 
will  not  launch  offensive  missiles  during  a  Soviet  attack 
and  will  ride  out  the  first  strike. 

It  is  assumed  that  the  Soviets  will  attack  the  MX  com¬ 
plex  uniformly  (i.e.,  an  equal  number  of  RVs  will  be  targeted 
at  each  complex  and  these  RVs  will  be  evenly  distributed 
among  the  shelters  in  each  complex). 

Methodology 

Modeling  Phase.  The  system  science  paradigm  was  used  to 
develop  the  models,  which  means  that  the  models  were  developed 
through  an  iterative  process  of  conceptualization,  analysis 
and  measurement,  and  computerization.  This  iterative  process 
was  an  aid  in  obtaining  all  the  desired  parameters  in  the 
models  (Ref  31) . 
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The  models  are  programmed  in  the  simulation  language 
Q-GERT  developed  by  A.  Alan  B.  Pritsker.  Q-GERT  was  chosen 
because  it  is  based  upon  queuing  theory,  and  the  entire 
system  being  modeled  is  envisioned  as  a  queuing  network  with 
the  attacking  RVs  and  defense  systems  representing  the  cus¬ 
tomers  and  servers,  respectively.  Q-GERT  also  allows  Fortran 
to  be  directly  inserted  and  it  readily  utilizes  probabilities 
(Ref  26). 

Analysis  Phase.  The  model  output  provides  data  to 
develop  an  equation  (using  regression  analysis)  in  order  to 
determine  the  following: 

1.  The  number  of  shelters  and  the  number  of  exo- 
atmospheric  interceptors  which  provide  a  prescribed 
level  of  survivability  with  one  terminal  defense 
unit  (TDU) . 

2.  The  number  of  shelters  which  provide  a  prescribed 
level  of  survivability  with  one  TDU. 

3.  The  number  of  shelters  which  provide  a  prescribed 
level  of  survivability  with  two  TDUs. 

With  these  data  and  equations  in  hand,  each  alternative  was 
quantitatively  (i.e.,  in  dollars)  examined,  and  conclusions 
were  drawn  based  upon  both  quantitative  and  qualitative 
(ex.,  political  and  environmental)  considerations.  An  out¬ 
growth  of  this  analysis  will  be  the  significant  variables 
of  this  study  and  recommendations  for  further  study  in  the 
area  of  Ballistic  Missile  Defense  (BMD) . 


II.  The  Model 


The  system  being  modeled  can  be  divided  into  four 
separate  subsystems:  attack,  target,  exoatmospheric  defense, 
and  endoatmospheric  defense  (Figure  2).  The  variables  and 
parameters  of  each  subsystem,  the  values  assigned  to  the 
variables,  and  sensitivity  analysis  on  these  values  will  be 
discussed  in  this  chapter;  as  well  as  the  probability  of 
kill  methodology,  the  simulation  models,  and  the  verification 
and  validation  of  the  models. 

Attack  Subsystem 

The  attack  subsystem  is  composed  of  Soviet  ICBMs  which 
transport  the  attacking  RVs  designed  to  destroy  a  hardened 
MX  shelter.  By  the  mid-1980s,  the  Soviet  Union  is  expected 
to  have  approximately  6000  warheads  in  the  one -megaton  yield 
range  (Refs  28:23  and  26).  It  is  believed  that  these  war¬ 
heads  could  be  delivered  in  various  configurations  by  the 
SS-17,  SS-18,  SS-19,  and  older  generation  Soviet  ICBMs.  The 
CEP  of  these  warheads  ranges  from  0.1  nautical  miles  (NM) 
to  0.2S  NM  (Refs  16:54;  28:23;  and  26).  By  decreasing  the 
yield  of  each  warhead,  thus  allowing  more  warheads  per 
missile,  the  number  of  attacking  RVs  could  be  increased  with¬ 
out  increasing  the  number  of  ICBMs.  This  process  of  fraction- 
ization  causes  fratracide  to  become  a  problem.  Fratracide  is 
the  destruction  of  an  RV  by  another  RV’s  detonation  (Ref  9:34). 
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ATTACK  SUBSYSTEM 


There  are  basically  two  kinds  of  fratracide,  local  and  area. 
Local  and  area  fratracide  deal  respectively  with  the  deton¬ 
ation  of  an  earlier  warhead  assigned  to  the  same  target  and 
the  detonation  of  an  earlier  warhead  assigned  to  a  different, 
though  nearby,  target . (Ref  2:5S).  In  a  fratracide  model 
developed  by  Steinbruner  and  Garwin,  the  second  warhead's 
survivability  level  varies  and  the  third  and  fourth  warheads 
usually  will  not  survive  (Ref  35).  In  their  model,  the 
authors  assume  that  individual  warheads  arriving  at  the  same 
target  must  be  separated  by  at  least  six  minutes  or  they  have 
a  high  probability  of  being  destroyed  by  the  first  detonation. 
This  thesis  uses  RV  spacing  of  from  several  seconds  to  slightly 
over  two  minutes.  "Normally  it  is  assumed  that  at  most  two, 
and  perhaps  only  one,  warhead  can  be  used  per  target  without 
being  overcome  by  fratracide"  (Ref  2:58). 

Target  Subsystem 

The  target  subsystem  consists  of  200  MX  missiles  located 
in  from  1600  to  4600  hardened  horizontal  MX  shelters,  and 
any  terminal  defense  units  (TDUs)  and  its  components  (Figure  3). 
The  one  MX  missile  hidden  in  an  MX  complex  will  represent  from 
8  to  23  aimpoints,  depending  on  the  number  of  shelters  in  the 
complex.  The  destruction  of  more  than  one  shelter  by  a  single 
attacking  RV  is  prevented  by  spacing  the  shelters  approximately 
5200  feet  apart  (Ref  15).  The  hardness  of  each  horizontal 
shelter  is  dependent  upon  the  shelter  design  selected.  With 
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four  roof  portals  (to  allow  Soviet  verification),  the  expected 
hardness  of  an  MX  or  TDU  shelter  to  overpressure  would  be 
approximately  600  pounds  per  square  inch  (psi)  (Ref  19).  A 
decrease  in  the  number  and/or  size  of  the  roof  portals,  however, 
could  increase  the  expected  hardness  to  over  1000  psi  (Ref  16:58) 

The  MX  and  TDU  shelters  are  also  vulnerable  to  the  crater¬ 
ing  caused  by  an  attacking  Soviet  RV.  Appendix  A  provides  a 
diagram  of  a  horizontal  shelter  and  the  methodology  for  calcu¬ 
lating  the  vulnerability  to  cratering  for  both  a  MX  and  TDU 
shelter . 

If  a  terminal  defense  system  exists,  its  radar  network 
becomes  part  of  the  target  subsystem.  The  TDU  will  not  be 
able  to  defend  the  MX  if  the  radar  network  is  destroyed,  but 
this  consequence  is  overcome  since  it  is  felt  that  the  TDU's 
radar  network  can  be  made  survivable  (Ref  30). 

Exoatmo spheric  Defense  Subsystem 

This  defensive  subsystem  is  currently  envisioned  as 
having  two  phases:  (1)  data  gather ing/attack  evaluation,  and 
(2)  intercept.  The  attack  data  will  come  from  a  ballistic 
information  probe  launched  immediately  when  alerted  to  a 
possible  Soviet  attack  by  other  national  resources.  At  an 
altitude  of  approximately  300,000  feet  the  probe  will  scan  the 
attack  azimuth  with  long  wave  infrared  (LWIR)  sensors  and, 
using  on-board  data  processing,  evaluate  the  attack,  sending 
the  information  back  to  command  authorities.  At  this  altitude, 
the  information  probe  can  view  the  approaching  vehicles 


against  a  cold  space  background  -  a  requirement  for  accurate, 
timely  IR  detection.  This  phase  is  beyond  the  problem  being 
studied  and  is  not  modeled  or  discussed  hereafter. 

The  intercept  phase  occurs  when  a  ballistic  intercept 
probe  carrying  multiple  high-speed  interceptors  is  launched. 

A  LWIR  sensor  on  the  intercept  probe  detects  individual  RVs 
and  assigns  each  interceptor  missile  it  is  carrying  to  an  RV. 
Each  interceptor  then  uses  LWIR  terminal  guidance  and  will 
directly  impact  the  incoming  RV  to  destroy  it.  This  kill 
mechanism  is  very  similar  to  that  planned  for  the  current 
F-15  launched  antisatellite  (ASAT)  program. 

"Impact -point  prediction  ...  is  not  expected  to  be 
able  to  resolve  impacts  among  closely  spaced  sh elters  of  any 
of  the  MPS  emplacement  schemes  under  coins  ider  at  ion"  (Ref  1:7). 
Therefore,  the  exoatmospheric  BMD  system  cannot  preferentially 
defend  the  MX  or  terminal  defense  unit,  and  must  intercept 
each  attacking  RV  that  comes  within  the  defense  systems'  field 
of  view.  The  assumption  that  all  interceptors  can  be  fired 
while  the  entire  attack  cloud  is  scanned,  coupled  with  the 
lack  of  preferential  defense  capability,  results  in  the  exo¬ 
atmospheric  defense  system  firing  its  N  interceptors  at  the 
first  N  attacking  RVs  or  a  one-to-one  basis.  Thus,  the 
defense  system  is  defending  itself  by  not  allowing  any  RV  to 
come  close  enough  to  destroy  it. 

The  circular  error  probable  (CEP)  of  the  exoatmospheric 
interceptors  is  assumed  to  be  normally  distributed  in  both 
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the  horizontal  and  vertical  dimensions  of  the  attacking  RV 
plane.  Since  it  is  assumed  that  the  exoatmospheric  inter¬ 
ceptor  will  pass  through  the  plane  perpendicular  to  the  RV 
trajectory,  the  third  dimension  of  depth  is  ignored  and 
circular  error  probable  rather  than  spherical  error  probable 
calculations  are  made. 

Endoatmospheric  Defense  Subsystem 

The  endoatmospheric  defense  subsystem  includes  the 
terminal  defense  unit's  (TDU's)  radar  network,  three  hypersonic 
nuclear  armed  interceptors,  and  a  control  unit.  The  control 
unit,  the  interceptors,  and  part  of  the  radar  network  will  be 
located  in  a  hardened,  horizontal  shelter  (Ref  12). 

The  radar  network  of  this  subsystem,  operating  without 
the  exoatmospheric  defense  layer,  might  have  three  stages. 

The  first  stage  would  be  an  early-warning  system  which  would 
detect  attacking  RVs  targeted  somewhere  within  the  200  MX 
missile  field.  This  stage  would  not  be  needed  if  the  endo¬ 
atmospheric  defense  system  were  operating  as  part  of  a  layered 
system.  The  second  stage  might  be  an  MX  complex  radar  warning 
system.  This  stage  would  detect  attacking  RVs  targeted  on  a 
particular  MX  complex.  If  attacking  RVs  are  descending  on  the 
complex,  the  last  stage  of  the  radar  network  would  begin  to 
function.  This  radar  would  track  incoming  RVs  and  determine 
their  precise  target  among  the  shelters  of  one  MX  complex 


■i  *«wii  wgi-Qta 


>  *•**•»*-«  Ho**-  !***A**  y< 


The  strategy  selected  for  the  terminal  defense  unit 
(TDU)  allows  a  TDU  to  launch  the  first  two  interceptors  at 
RVs  aimed  at  either  the  MX  or  TDU  shelters.  The  remaining 
interceptor  will  be  used  for  MX  defense  only.  Therefore,  if 
an  RV  is  attacking  a  TDU  shelter  and  only  one  interceptor 
remains,  the  interceptor  will  not  be  launched,  and  the  TDU 
will  be  subjected  to  an  RV  detonation  which  may  or  may  not 
destroy  the  TDU  (Ref  22:38).  This  strategy  was  selected 
because  it  was  determined  most  effective  in  James  Moore’s 
recent  study  (Ref  22). 

The  control  unit  of  the  TDU  determines  whether  or  not 
an  interceptor  should  be  launched  based  upon  the  above 
strategy.  The  process  of  launching  an  interceptor  requires 
the  TDU  to  leave  the  shelter,  acquire  an  attacking  RV  with 
its  radar,  launch  an  interceptor,  and  return  to  the  shelter 
(Ref  22:13).  All  endoatmospheric  interceptors  are  assumed 
to  detonate  at  approximately  20,000  feet  and  attempt  to 
destroy  an  attacking  RV  with  neutron  radiation.  Circular 
error  probable  (CEP) ,  rather  than  spherical  error  probable 
(SEP),  is  used  since  it  is  assumed  that  the  interceptor 
passes  through  the  RV  plane  as  in  the  exoatmospheric  defense 
layer . 

System  Variables 

The  subsystems  previously  described  contain  the  following 
set  of  variables  which  have  a  direct  or  indirect  effect  on 
MX  survivability.  * 


t 


17 


The  set  of  attack  subsystem  variables  which  can  affect 
MX  survivability  are  the  number  of  RVs  attacking  an  MX  com¬ 
plex,  the  Soviet  targeting  strategy,  and  the  yield  and  CEP 
of  the  attacking  RVs.  Since  the  RVs  may  be  subjected  to  a 
shower  of  neutrons,  their  ability  to  survive  the  neutron 
fluence  of  the  interceptor  warhead  is  a  major  factor  in 
determining  MX  survivability.  The  height  of  burst  affects 
an  attacking  RV's  overpressure  and  cratering  capability  which 
also  have  a  direct  impact  on  MX  survivability.  RV  reliabil¬ 
ity  is  also  an  important  factor. 

There  are  various  target  subsystem  variables  which  can 
affect  MX  survivability.  Shelter  hardness,  which  depends  on 
shelter  design  and  soil  type,  is  a  variable  considered  to 
have  a  definite  effect  on  MX  survivability.  The  number  of 
shelters  per  complex,  shelter  spacing,  and  target  altitude  are 
other  target  subsystem  variables  which  must  be  considered 
when  determining  MX  survivability. 

Exoatmospheric  defense  variables  are  the  weapon  radius 
(WR)  and  CEP  of  the  interceptors,  the  number  of  interceptors, 
and  the  reliability  of  the  infrared  guidance  system. 

The  set  of  endoatmospheric  defense  variables  which  can 
impact  the  survivability  of  the  MX  include  the  number  of 
TDUs  per  complex,  the  number  of  interceptors  per  TDU,  and 
the  launch  strategy  of  the  TDU.  Other  important  variables 
include  the  yield,  CEP,  and  reliability  of  the  interceptors; 
and  the  reliability  of  the  TDU’s  radar  network. 
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The  interaction  of  all  of  the  above  variables  determines 
the  survivability  of  U.S.  MX  missiles.  The  residual  effects 
of  a  nuclear  detonation  (ex.,  fallout)  and  the  treaty  banning 
above  ground  nuclear  testing  preclude  using  actual  nuclear 
weapons  to  study  MX  survivability.  Therefore,  this  study 
will  use  models  which  employ  all  variables  considered  import¬ 
ant  in  determining  the  outcome  of  a  Soviet  attack  on  U.S. 

MX  complexes. 

Structural  Model 

The  variables  previously  mentioned  were  selected  because 
they  were  deemed  to  have  a  significant  effect  on  MX  surviva¬ 
bility.  Estimates  of  MX  survivability  are  critical  to  attain¬ 
ing  the  objectives  of  this  thesis.  Some  of  these  variables 
are  given  preassigned  values,  while  others  will  be  parameters 
in  the  models  and  will  be  varied  on  each  run  of  the  model. 

The  variables  chosen  for  modeling  the  attack  subsystem 
are  the  number  of  Soviet  RVs  attacking  an  MX  complex;  the 
targeting  strategy;  the  yield,  CEP,  and  height  of  burst  of 
the  RV;  RV  reliability;  and  the  sure-safe  and  sure-kill  neutron 
fluence  of  the  attacking  RV.  Of  these  variables,  only  the 
number  of  attacking  RVs  is  considered  a  parameter. 

The  target  subsystem  variables  chosen  for  the  modeling 
phase  of  this  study  are  the  number  of  shelters  per  complex, 
the  sure-safe  and  sure-kill  overpressure  levels  of  a  shelter, 
the  type  of  soil  in  which  the  shelters  are  constructed,  and 
the  altitude  of  the  area  in  which  the  shelters  are  located. 
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The  number  of  shelters  per  complex  is  the  only  model 
parameter. 

The  variables  included  in  modeling  the  exoatmospheric 
defense  subsystem  are  the  number  of  interceptors;  the 
strategy,  weapon  radius,  and  CEP  of  the  interceptors;  and 
the  reliability  of  the  infrared  (IR)  guidance  system.  The 
number  of  exoatmospheric  interceptors  is  the  only  exoatmos¬ 
pheric  defense  model  parameter. 

The  variables  chosen  for  modeling  the  endoatmospheric 
defense  subsystem  are  the  number  of  TDUs  per  complex;  the 
number  of  interceptors  per  TDU;  interceptor  strategy,  yield, 
CEP,  and  reliability;  and  the  radar  network  reliability.  The 
number  of  TDUs  per  complex  is  the  only  parameter  selected 
for  the  modeling  of  the  endoatmospheric  defense  system. 

The  variables  which  are  not  considered  parameters  and 
are  given  preassigned  values  establish  limitations  on  the 
results  of  this  thesis.  The  assumed  targeting  strategy  is 
to  randomly  target  each  shelter  until  all  shelters  are  tar¬ 
geted,  and  then  randomly  assign  excess  RVs  to  those  same 
shelters.  If  there  are  twice  as  many  RVs  as  shelters,  then  the 
Soviets  will  randomly  target  each  shelter  on  a  two-to-one  basis 
If  there  are  fewer  RVs  than  shelters,  the  RVs  are  just  ran¬ 
domly  targeted  against  the  shelters.  The  two  main  ways  to 
destroy  a  ground  emplaced  shelter  are  with  overpressure,  or 
cratering  and  ground  shock.  The  optimum  height  of  burst 
(HOB)  for  overpressure  depends  on  weapon  yield  and  could  be 
several  thousand  feet  above  ground  level.  The  optimum  HOB 


for  cratering  depends  on  weapon  yield  and  soil  type  and  could 
be  several  hundred  feet  below  ground  level.  In  the  models, 
the  attacking  RVs  detonate  when  they  contact  the  earth's 
surface  (i.e.,  height  of  burst  equal  to  zero).  While  this 
height  of  burst  is  not  optimum  for  either  overpressure  or 
cratering,  it  allows  for  sizeable  contributions  to  shelter 
destruction  from  each  effect  and  facilitates  evaluation  of 
both  kill  mechanisms.  It  also  allows  a  two-dimensional  (CEP) 
analysis  instead  of  a  three-dimensional  (SEP)  analysis  of  the 
RV-shelter  interaction.  The  yield  of  the  RVs  will  be  one 
megaton  since  the  Soviets  will  have  approximately  6000  war¬ 
heads  in  this  range  by  1985  (Ref  28:23  and  26).  The  CEP  of 
the  Soviet  SS-18  is  within  the  rang®  of  0.12  to  0.25  NM 
(Ref  16:54),  and  for  the  purpose  of  this  study  is  assumed  to 
be  0.2  NM.  Reasonable  estimates  for  the  sure-kill  and  sure- 
safe  neutron  fluence  of  the  RVs  have  been  established  at  1017 
and  1013  neutrons  per  square  centimeter  (N/CM2),  respectively 
(Ref  7).  The  reliability  of  the  incoming  RVs  has  been  assumed 
to  be  one. 

Approximate  levels  for  the  sure-kill  and  sure-safe  over¬ 
pressure  of  the  MX  and  TDU  shelters  have  been  set  at  1250 
and  750  pounds  per  square  inch  (psi),  respectively  (Ref  22:36). 
The  soil  type  in  which  the  MX  and  TDU  shelters  are  built  is 
assumed  to  be  dry  soil  and/or  dry  soft  rock,  representative 
of  the  West  or  Southwest  U.S.  where  the  shelters  might  be 
constructed.  The  shelters  are  assumed  to  be  spaced  suffici- 
ciently  far  apart  to  prevent  one  Soviet  RV  from  destroying  two 
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shelters.  Each  of  the  shelters  is  assumed  to  be  at  an 
altitude  of  4000  feet  (MSL)  ,  which  is  a  representative 
altitude  for  the  valleys  of  Utah  and  Nevada. 

The  diameter  of  the  exoatmospheric  interceptor  is 
equal  to  that  of  the  current  antisatellite  (ASAT)  design 
(i.e.,  one  foot)  (Ref  10:244).  This  one  foot  diameter, 
coupled  with  a  one  foot  wide  mesh  which  will  unfold  prior 
to  impact,  results  in  a  one  and  one -half  foot  weapons  radius 
(Figure  4).  Since  it  is  assumed  that  the  attacking  RVs  are 
cone  shaped  with  a  radius  of  approximately  one  and  one -half 
feet,  the  effective  weapons  radius  of  the  exoatmospheric 
interceptor  is  set  at  three  feet  (Figure  4).  Since  an 
unclassified  value  for  long  wave  IR  guidance  accuracy  was 
unavailable,  the  CEP  of  the  exoatmospheric  interceptor  was 
chosen  to  be  two  feet  based  upon  the  resolution  accuracy  of 
a  Department  of  Defense  (DoD)  infrared  telescope  currently 
in  development  (Ref  29:24),  and  the  general  capabilities  of 
current  air-to-air  missiles.  The  strategy  of  the  exoatmos¬ 
pheric  defense  system  is  to  launch  all  of  its  N  interceptors 
at  the  first  N  RVs  attacking  the  U.S.  MX  complex  on  a  one- 
to-one  basis.  The  reliability  of  the  interceptors  and  the  IR 
guidance  system  have  been  assumed  to  be  one. 

The  number  of  endoatmospheric  interceptors  per  TDU  has 
been  taken  as  three  based  on  current  DoD  plans  (Ref  37). 

The  yield  of  the  nuclear-armed  endoatmospheric  interceptors 
is  set  at  five  kilotons  (KTs),  since  it  is  assumed  that  the 
U.S.  prefers  to  detonate  as  small  a  nuclear  warhead  as  possible 
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Figure  4.  Exoatmospheric  Interceptor  and  RV 

within  its  atmosphere.  A  circular  error  probable  (CEP) 
of  600  feet  is  selected  as  being  representative  of  U.S. 
technological  capabilities  (Ref  22:21).  This  CEP  is  much 
larger  than  the  exoatmospheric  interceptor's  because  the 
endoatmospheric  interceptor  does  not  use  terminal  guidance 
and  must  rely  on  more  conventional  guidance  techniques. 

Long  wave  IR  terminal  guidance  is  currently  not  feasible 
in  the  endoatmospheric  intercept  phase  due  to  the  background 
interference  generated  by  the  atmosphere  at  these  wavelengths 
(approximately  10’^  meters).  As  mentioned  previously,  the 
strategy  of  the  TDU  allows  launch  of  the  first  two  interceptors 
at  RVs  aimed  at  either  the  MX  or  TDU  shelters,  and  the  remain- 
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ing  interceptor  to  be  used  for  MX  defense  only  (Ref  22:38). 

The  interceptors  and  the  radar  network  are  assumed  to  be 
100%  reliable. 

The  ballistic  missile  defense  models  allow  the  number 
of  shelters  per  complex,  the  number  of  attacking  RVs,  the 
number  of  exoatmospheric  interceptors,  and  the  number  of 
TDUs  per  complex  to  be  selected  in  order  to  determine  a  level 
of  MX  survivability. 

Sensitivity  of  Model  Variables 

A  few  of  the  model  variables  which  are  given  preassigned 
values  are  very  critical  in  determining  the  effectiveness  of 
the  exoatmospheric  interceptors,  the  endoatmospheric  inter¬ 
ceptors,  and  the  Soviet  RVs;  all  of  which  have  an  impact  on 
MX  survivability.  Precise  values  of  these  variables  are 
classified  or  unknown,  which  necessitates  sensitivity  analysis 

The  critical  variables  which  determine  exoatmospheric 
interceptor  probability  of  kill  (PK)  are  the  interceptor  CEP 
and  the  effective  weapons  radius  (interceptor  radius  plus  RV 
radius)  of  the  interceptor.  The  PKs  of  individual  inter¬ 
ceptors  for  various  CEP  and  effective  weapons  radius  (WR) 
combinations  were  calculated  as  shown  in  Appendix  E  and  are 
presented  in  Table  I,  Since  the  radius  of  Soviet  RVs  is 
beyond  U.S.  control,  the  changes  in  effective  weapon  radius 
are  essentially  changes  in  interceptor  weapon  radius.  An 
effective  weapon  radius  larger  than  3.5  feet  was  not  examined 


TABLE  I 

Exoatmospher ic  Interceptor 

PKs 

Effective 

WR  (feet) 

0.5 

1.0 

CEP  (feet) 

l.S  2.0 

2.5 

3.0 

2.0 

1.0 

.94 

.71  .50 

.36 

.27 

2.5 

1.0 

.99 

.85  .66 

.50 

.38 

3.0 

1.0 

1.0 

.94  .79 

.63 

.50 

3.5 

1.0 

1.0 

.98  .88 

.74 

.61 

because  it  was  felt  that  interceptors  of  this  size/mass 
would  not  be  able  to  maneuver  fast  enough  at  the  high 
speeds  experienced  during  an  RV  engagement.  In  general, 
a  one-half  foot  decrease  in  CEP  provides  a  greater  increase 
in  interceptor  PK  than  a  one-half  foot  increase  in  inter¬ 
ceptor  weapon  radius.  Therefore,  interceptor  CEP  is  considered 
more  critical  than  weapon  radius. 

Interceptor  CEP  and  yield  are  the  interceptor  variables 
which  are  critical  in  determining  the  PK  of  endoatmo spheric 
interceptors.  Various  combinations  of  interceptor  CEP  and 
yield,  and  their  corresponding  PKs  are  shown  in  Table  II.  The 
PKs  were  calculated  using  the  neutron  kill  subroutines  of  the 
simulation  program  (Appendix  G).  This  table  shows  that 
interceptor  PK  increases  at  a  decreasing  rate  when  yield  is 
increased  and  CEP  is  held  constant,  and  interceptor  PK 
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TABLE  II 


Endoatmospheric  Interceptor  PKs 
(1013  N/CM2  Sure-Safe,  1017  N/CM2  Sure-Kill) 


Yield 

(KT) 

200 

CEP 

400 

(feet) 

600 

800 

2.5 

.75 

.60 

.49 

.42 

5.0 

.83 

.70 

.60 

.52 

10.0 

.89 

.79 

.70 

.62 

15.0 

.92 

.84 

.76 

.68 

20.0 

.94 

.86 

.79 

.71 

increases  at  a  fairly  constant  rate  when  CEP  is  decreased 
and  yield  is  held  constant.  Therefore,  as  in  the  exoatmos- 
pheric  interceptor  case,  CEP  is  considered  the  most  critical 
interceptor  variable  in  determining  interceptor  PK.  Increas¬ 
ing  the  yield  of  an  interceptor  may  also  be  unacceptable 
because  larger  nuclear  warheads  detonated  in  the  atmosphere 
cause  increased  nuclear  fallout  and  peripheral  blast  damage. 

Other  important  variables  which  impact  the  endoatmos¬ 
pheric  interceptors'  PK  are  the  sure-safe  and  sure-kill 
neutron  fluence  levels  of  the  attacking  RVs.  Various  combin¬ 
ations  of  sure-safe  and  sure-kill  fluence  levels  are  presented 
in  Table  III  for  the  interceptor  yield  and  CEP  used  in  this 
thesis.  The  PKs  in  Table  III  were  calculated  using  the  neutron 
kill  subroutines  of  the  simulation  program  (Appendix  G).  Since 
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the  models  in  this  study  use  sure-safe  and  sure-kill  neutron 
intensity  levels  of  1013  and  1017  N/CM2,  respectively,  this 
table  shows  that  varying  either  the  sure-safe  or  sure-kill 
fluence  has  a  significant  effect  on  interceptor  PK. 


TABLE  III 

Endoatmospheric  Interceptor  PKs 
(CEP  »  600  feet,  Yield  *  5  kilotons) 


Sure -Kill 
Fluence  (N/CM2) 

Sure- 

101 1 

Safe  Fluence  (N/CM2) 

1012  1013  1014 

1015 

.98  5 

.  984 

.983 

.973 

1018 

.915 

.882 

.813 

.634 

1017 

.800 

.726 

.599 

.390 

1018 

.680 

.583 

.447 

.277 

The  critical  variables  that  determine  an  attacking  RV's 
PK  are  the  yield  and  CEP  of  the  RV,  and  the  sure- safe  and 
sure-kill  overpressure  levels  of  the  MX  or  TDU  shelter  being 
attacked.  The  PKs  of  a  single  RV  are  shown  in  Tables  IV, 

V,  and  VI  for  various  combinations  of  RV  CEP  and  yield,  and 
MX  shelter  sure- safe  and  sure-kill  overpressure  levels.  These 
PKs  were  calculated  using  the  overpressure  kill  routines  in 
the  simulation  program  (Appendix  G) .  In  general,  these  tables 
show  that,  for  yields  greater  than  about  1000  kilotons,  an 
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TABLE  IV 

PK  Against  MX  Shelter 


( 500  psi  Sure-Safe,  1000  psi  Sure-Kill) 


Yield 

(KT) 

0.10 

CEP 

0.14 

(NM) 

0.20 

0.25 

250 

.965 

.599 

.549 

.190 

500 

1.0 

.751 

.  572 

.543 

750 

1.0 

.986 

.586 

.557 

1000 

1.0 

1.0 

.601 

.566 

1250 

1.0 

1.0 

.648 

.574 

1500 

1.0 

1.0 

.763 

.581 

Yield 

(KT) 


250 

500 

750 

1000 

1250 

1500 


TABLE  V 

PK  Against  MX  Shelter 


(750  psi  Sure-Safe,  1250  psi  Sure-Kill 


0.10 


.739 

1.0 

1.0 

1.0 

1.0 

1.0 


CEP  (NM) 

0.14  0.20 


.599 

.634 

.805 

.992 

1.0 

1.0 


.442 

.572 

.586 

.599 

.609 

.619 


0.25 


.160 

,455 

.557 

.566 

.574 

.581 
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TABLE  VI 

PK  Against  MX  Shelter 
(900  psi  Sure-Safe,  1400  psi  Sure-Kill) 


Yield 

(KT) 

0.10 

CEP 

0.14 

(NM) 

0.20 

0.25 

250 

.677 

.599 

.271 

.161 

500 

1.0 

.635 

.572 

.  282 

750 

1.0 

.678 

.  586 

.  548 

1000 

1.0 

.940 

.599 

.566 

1250 

1.0 

.999 

.609 

.574 

1500 

1.0 

1.0 

.619 

.  581 

increase  in  yield  does  not  have  a  significant  effect  on  RV 
PK;  and  the  sure-safe  and  sure-kill  overpressure  levels  are 
not  critical  in  determining  RV  probability  of  kill.  The 
most  significant  effect  on  RV  PK  occurs  when  the  CEP  of  the 
reentry  vehicles  change  from  0.2  NM  to  0.14  NM,  and  vice 
versa.  Considering  that  the  models  use  sure-safe  and  sure- 
kill  overpressure  levels  of  750  psi  and  1250  psi,  respectively, 
and  a  CEP  and  yield  of  0.2  NM  and  1000  kilotons,  respectively, 
the  attacking  RVs  CEP  is  considered  the  most  critical  variable 
in  determining  the  probability  of  killing  the  MX. 
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Probabilities  of  Kill 

In  order  to  perform  the  probability  of  kill  (PK) 
calculations  for  the  exoatmospheric  interceptors,  the  TDU 
interceptors,  and  the  Soviet  RVs,  various  assumptions  must 
be  made.  It  is  assumed  that  the  exact  MX  and  TDU  locations 
within  the  MX  complex  are  unknown  and  thus  the  attacking  RVs 
are  randomly  targeted  among  the  entire  complex  of  shelters. 

This  and  the  fact  that  the  exoatmospheric  defense  system 
cannot  preferentially  defend  the  MX  complex  make  exoatmos¬ 
pheric  intercept  strategy  immaterial  to  the  model  results. 
Another  assumption  is  that  the  exoatmospheric  defense  system 
has  the  time  and  capability  to  scan  the  entire  RV  attack 
cloud  while  all  the  interceptors  are  being  launched.  The 
exoatmospheric  interceptors  have  a  direct  impact  kill  mechan¬ 
ism  with  an  effective  "cookie-cutter"  weapon  radius  (WR)  of 
three  feet. 

The  products  of  the  detonation  of  nuclear-armed  endo- 
atmospheric  interceptors  in  the  atmosphere  are  overpressure, 
dynamic  pressure,  thermal  radiation,  gamma  rays,  and  neutron?. 

A  reentry  vehicle  (RV)  is  designed  to  withstand  the  high 
temperatures  which  occur  when  it  reenters  the  earth's  atmos¬ 
phere,  hence  it  is  probably  capable  of  withstanding  the  effects 
of  thermal  radiation.  The  effects  of  overpressure  and  dynamic 
pressure  also  have  little  impact  on  attacking  RVs  because  of 
their  aerodynamic  "low  drag"  design  (Ref  6).  Gamma  rays  are 
not  an  effective  kill  mechanism  because  the  prompt  (source) 
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gammas  are  such  a  small  percentage  of  the  total  energy  of 
the  explosion  that  the  intensity  levels  are  not  sufficiently 
high  to  cause  the  heating  required  to  crack  or  disrupt  the 
fissile  material  in  the  RV.  Secondary  gamma  fluence  is  also 
too  low  because  secondary  gammas  are  produced  by  neutron-air 
reactions  over  the  entire  volume  of  air  populated  by  the 
neutrons  and  not  from  a  point  source  (the  weapon)  (Ref  4). 

The  neutrons  created  by  the  detonation  of  an  interceptor, 
however,  can  destroy  an  RV  by  heating  the  fissile  material 
if  the  neutron  fluence  is  sufficiently  high  (Ref  11:1137). 
Therefore,  the  endoatmospheric  interceptors'  kill  mechanism 
is  the  neutrons  created  by  their  detonation.  Appendix  D 
presents  the  method  for  calculating  endoatmospheric  proba¬ 
bility  of  kill  (PK)  based  upon  neutron  fluence. 

The  effects  of  a  surface  burst  nuclear  explosion  are 
dynamic  pressure,  overpressure,  neutrons,  gamma  rays,  thermal 
radiation,  ground  motion,  and  cratering.  The  thick,  steel 
and  reinforced  concrete  walls  of  a  shelter  provide  shielding 
against  gamma  rays,  neutrons,  and  thermal  radiation  (Ref  18). 
The  MX  shelters  will  also  be  built  on  a  suspension  system  to 
prevent  damage  from  ground  motion  (Ref  33).  If  a  shelter  is 
built  flush  with  the  ground,  the  destructive  sideloadings  of 
dynamic  pressure  can  also  be  avoided  (Ref  22:17).  Although 
a  shelter  can  be  designed  to  limit  damage  by  overpressure, 
a  sufficiently  high  level  of  overpressure  will  destroy  the 
shelter.  A  shelter  can  also  be  destroyed  if  the  shelter  is 
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within  the  crater  radius  caused  by  the  nuclear  explosion, 
or  rendered  inoperable  if  the  ejecta  caused  by  the  crater 
is  thick  enough  to  cover  the  door  of  a  horizontal  shelter. 
Hence,  a  horizontal  shelter  can  be  destroyed  or  rendered 
inoperable  by  the  effects  of  overpressure  and  cratering. 
Appendices  A  and  C  show  the  model  procedures  for  computing 
the  Soviet  RV  probability  of  kill  due  to  cratering  and  over¬ 
pressure,  respectively. 

The  Simulation 

The  models  in  this  thesis  simulate  the  defense  of  an  MX 
complex  being  attacked  by  Soviet  RVs  using  the  simulation 
language  Q-GERT.  The  computer  programs  and  Q-GERT  networks 
are  shown  in  Appendix  G.  The  BUS  deployment  and  defense 
service  times  were  included  to  provide  realistic  timing 
control  of  the  model  flow,  but  they  do  not  have  a  significant 
effect  on  the  model  results.  It  is  assumed  that  the  defense 
systems  cannot  be  saturated. 

The  simulation  first  generates  the  required  number  of 
attacking  RVs  which  are  dispensed  from  the  RV  BUS  normally 
distributed  in  time  with  a  mean  of  five  seconds  and  a 
standard  deviation  of  two  seconds.  If  the  number  of  RVs  is 
exactly  twice  the  number  of  shelters,  two  RVs  are  randomly 
assigned  to  each  shelter.  The  number  of  RVs  will  never  be 
greater  than  twice  the  number  of  shelters  due  to  the  fratri¬ 
cide  limit  previously  mentioned.  If  the  number  of  RVs  is 
larger  than  the  number  of  shelters  but  smaller  than  twice 
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the  number  of  shelters,  the  simulation  insures  that  a  minimum 
of  one  RV  is  targeted  at  each  shelter  if  the  Soviets  possess 
enough  RVs  to  target  each  shelter  at  least  once.  Excess  RVs 
are  also  randomly  targeted  on  the  shelters.  After  assigning 
all  the  attacking  RVs  to  shelters,  all  the  RVs  are  given  an 
eleven  minute  flight  time  delay  until  they  reach  the  exoatmos- 
pheric  defense,  if  one  exists. 

Layered  Defense  CExoatmospheric  Layer).  The  first  layer, 
the  exoatmospheric  defense,  attempts  to  defend  the  MX  complex 
by  launching  N  interceptors  on  a  one-to-one  basis  at  the  first 
N  RVs  encountered.  Although  the  RVs  are  intercepted  first-come 
first -serve,  the  exoatmospheric  defense  does  not  know  the 
intended  target  since  the  RVs  were  randomly  assigned  to  t  .e 
targets.  All  RVs  which  do  not  encounter  an  interceptor  or  are 
missed  by  an  interceptor  penetrate  the  exoatmospheric  defense 
layer  and  proceed  to  the  endoatmospheric  defense  layer  which 
consists  of  one  terminal  defense  unit  (TDU). 

Endoatmospheric  Layer  COne  TDU).  This  layer  attempts  to 
defend  the  MX  and  TDU  shelters  based  upon  the  designated 
interceptor  strategy.  If  the  TDU  is  destroyed,  the  MX  cannot 
be  defended  by  this  layer  of  defense.  The  attacking  RVs  are 
defended  against  one  at  a  time  with  a  service  time  of  five 
seconds.  The  simulation  computes  the  number  of  RVs  which 
impact  the  MX  shelter  and  calculates  the  probability  of  the 
Soviet  attack  killing  the  MX  using  the  following  equation: 
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where  PK  ■  probability  of  kill  for  the  entire  attack 

PK^  *  probability  of  kill  for  one  RV 
n  ■  number  of  RVs  impacting  the  MX  shelter. 

The  simulation  then  compares  this  PK  to  a  random  number  to 
determine  if  the  MX  is  destroyed. 

Endoatmospheric  Defense  (Two  TDUs).  This  simulation 
attempts  to  defend  the  MX  and  TDU  shelters  according  to  the 
following  strategy.  The  strategy  designates  one  of  the  TDUs 
as  the  primary  and  the  remaining  TDU  as  the  backup.  The 
primary  TDU  defends  the  MX  complex  until  it  is  destroyed  or 
has  launched  two  of  its  three  interceptors.  The  primary  TDU 
will  launch  its  first  two  interceptors  at  RVs  aimed  at  either 
the  MX  or  a  TDU  shelter,  but  will  save  its  remaining  inter¬ 
ceptor  for  defending  the  MX.  If  the  primary  TDU  is  destroyed 
or  has  launched  two  of  its  three  interceptors,  the  backup 
TDU  assumes  the  defense  role  if  it  has  not  been  destroyed.  If 
both  TDUs  are  destroyed,  the  MX  cannot  be  defended  by  this 
layer  of  defense.  If  the  backup  TDU  is  destroyed  or  has 
launched  all  three  of  its  interceptors,  the  primary  TDU 
reassumes  the  defense  role  if  it  has  not  been  destroyed. 

Both  TDUs  will  not  defend  a  TDU  shelter  which  has  been  des¬ 
troyed.  Either  TDU  defends  against  attacking  RVs  one  at  a 
time  with  a  service  time  of  five  seconds.  The  simulation 


computes  the  number  of  RVs  which  impact  the  MX  shelter  and 
calculates  the  probability  of  the  Soviet  attack  killing  the 
MX  using  the  same  equation  as  the  one  TDU  model.  The  simu¬ 
lation  then  compares  the  probability  of  kill  of  the  equation 
to  a  random  number  to  determine  whether  or  not  the  MX  is 
destroyed. 

Verification 

Although  the  exoatmospheric  and  endoatmospheric  defense 
models  are  combined  to  form  a  layered  defense  model,  they  are 
verified  separately  since  they  are  considered  independent  of 
one  another.  Therefore,  if  both  the  exoatmospheric  and  endo¬ 
atmospheric  models  are  verified,  the  layered  defense  model 
will  also  be  verified. 

All  defensive  models  were  verified  by  simulating  a  24  RV 
Soviet  attack  on  15  MX  shelters  and  comparing  the  model 
results  with  results  derived  analytically.  The  output  of 
one  model  run  is  MX  destroyed  or  not  destroyed.  This  type 
of  output  is  a  Bernoulli  trial  and  the  results  of  multiple 
Bernoulli  trials  can  be  characterized  by  the  binomial  distri¬ 
bution  (Ref  32:191).  A  binomial  distribution  can  be  approxi¬ 
mated  by  a  normal  distribution  if  the  number  of  runs  or 
trials  (n)  is  sufficiently  large,  and  the  probability  of  the 
MX  being  destroyed,  p,  is  close  to  one-half.  In  general, 
this  approximation  is  good  if  np  >  3  when  p  s  0.5  or 
n(l-p)  >  3  when  p  >  0.5  (Ref  39).  This  is  true  of  the 
models'  outputs  when  n  ■  1200  .  The  output  of  1200  simulation 


runs  is  p  ,  the  probability  that  the  MX  is  destroyed.  Hence, 
since  we  are  using  p  as  an  estimate  of  the  true  probability 
of  kill,  p  ,  we  can  be  (l-a)100t  confident  that  the  error 
of  the  estimate  will  be  less  than  a  specified  amount  e  (in 
decimal)  when  the  sample  size  is: 

Za/2  ? 


where  ^a/2  ^e  two-tailed  standardized  normal  statistic 
(Ref  38:212). 

In  order  to  test  the  null  hypothesis  that  the  PK  of  the 
MX  calculated  by  the  model  equals  the  analytical  PK,  a 
hypothesis  test  of  the  proportions  using  a  normal  distri¬ 
bution  is  used  since  the  sample  size  (n  *  1200)  is  suf¬ 
ficiently  large.  The  hypotheses  are: 

Hn:  Pn  *  P 


Hl:  P0  »*  P 


The  test  is: 


P  ’  z. 


<  P  <  P  ♦  z 


and  if  pQ  falls  within  this  acceptance  region,  then  it  is 

A 

assumed  that  pQ  *  p  ;  otherwise,  HQ  is  rejected  and  it  is 

A 

concluded  that  p  f  p  where 
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PK  determined  analytically 


P  * 

A 

p  ■  PK  determined  by  the  model 

Z  *  two-tailed  standardized  normal  statistic 
a/  i 

n  »  number  of  simulation  runs  ( Ref  38:209-213). 

Exoatmospheric  Defense  Model.  This  model  was  verified 
by  simulating  12  exoatmospheric  interceptors  against  a  24  RV 
attack.  The  yield  and  CEP  of  the  Soviet  RVs  were  assumed  to 
be  sufficient  to  provide  each  RV  with  a  100%  PK  if  not  des¬ 
troyed  by  an  interceptor.  The  analytical  calculation  of  an 
RV's  probability  of  killing  the  MX,  P0  »  was  performed  as 
outlined  in  Appendix  E,  and  results  in  a  PK  of  74.8%.  The 
1200  model  simulations  resulted  in  a  PK  of  74.9%,  which  is 
within  2.5%  of  the  true  mean.  The  95%  confidence  interval 
for  these  data  is  72.4%  <  p  <  77.4%  .  Since  the  analytical 
PK  falls  within  this  confidence  interval,  the  null  hypothesis 
cannot  be  rejected  and  the  exoatmospheric  defense  model 
functions  properly. 

Endoatmospheric  Defense  Model  (One  TDU) .  The  verification 
of  this  model  was  accomplished  using  1200  simulations  of  the 
Soviet  attack  on  the  MX  complex.  The  yield  and  CEP  of  the 
Soviet  RVs  are  set  at  one-megaton  and  0.2  nautical  miles, 
respectively.  The  analytical  probability  of  an  RV  killing 
the  MX,  pQ  ,  was  calculated  as  shown  in  Appendix  F,  and 
results  in  a  PK  of  44.4%.  The  1200  model  simulations  resulted 
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in  a  PK  of  42.21,  which  is  within  2.8%  of  the  true  mean. 

The  95%  confidence  interval  for  these  data  is  39.8%  <  p  <  45.4%. 
Hence,  since  the  analytical  PK  falls  within  this  confidence 
interval,  the  null  hypothesis  cannot  be  rejected  and  the  endo- 
atmospheric  defense  model  with  one  TDU  functions  properly. 

Layered  Defense  Model.  As  stated  previously,  the  exo- 
atmospheric  and  endoatmospheric  (one  TDU)  models  are  considered 
independent  of  one  another.  Therefore,  since  both  of  these 
models  have  been  verified,  the  layered  defense  model  has  been 
verified  and  functions  properly. 

Endoatmospheric  Defense  Model  (Two  TDUs).  The  only 
difference  between  this  model  and  the  one  TDU  model  are  the 
situations  in  which  the  primary  TDU  transfers  the  responsi¬ 
bility  of  defense  to  the  backup  TDU,  and  vice  versa.  Conse¬ 
quently,  only  the  situations  which  cause  this  transfer  of 
responsibility  need  to  be  verified.  The  verification  of 
these  transfers  was  accomplished  by  tracing  25  simulations 
of  the  Q-GERT  model.  Since  all  of  these  simulations  per¬ 
formed  the  transfer  of  defense  responsibility  successfully 
and  represent  all  possible  situations  that  might  occur,  the 
endoatmospheric  model  with  two  TDUs  functions  properly. 

Validation 

The  exoatmospheric  and  endoatmospheric  defense  systems 
that  these  models  portray  will  not  be  feasible  until  the 
mid  and  late  1980s,  respectively.  It  is  impossible  to  compare 


38 


the  behavior  of  each  model  with  the  behavior  of  an  actual 
system  which  has  not  been  built.  An  attempt  has  been  made 
to  include  those  variables  of  the  actual  systems  which  are 
anticipated  to  have  a  significant  effect  on  MX  survivability. 
Reasonable  assumptions  were  used  to  model  those  areas  where 
unclassified  data  were  not  available.  These  assumptions 
are  based  upon  an  extensive  search  of  unclassified  litera¬ 
ture  and  discussions  with  qualified  individuals  of  the  Aero¬ 
nautical  Systems  and  Armament  Divisions  of  Air  Force  Systems 
Command,  and  the  Air  Force  Institute  of  Technology.  All  of 
the  models  in  this  study  follow  directly  from  the  systems 
and  data  of  current  unclassified  sources,  and  within  the 
stated  limitations,  the  models  are  valid. 


III.  The  Analysis 


Experimental  Design 

The  models  developed  in  this  study  are  designed  to 
furnish  estimates  of  MX  survivability  when  an  MX  complex  is 
defended  by  a  layered  defense,  one  terminal  defense  unit 
(TDU) ,  or  two  TDUs.  With  these  models,  the  effects  of  exo- 
atmospheric  interceptors,  protective  shelters,  and  TDUs  on 
MX  survivability  can  be  examined. 

Layered  Defense  Model.  The  experimental  design  of  this 
model  consists  of  two  factors  (number  of  shelters  and  number 
of  interceptors)  and  one  response  (MX  survivability)  as  shown 
in  Figure  5. 


HUMBER 
OF 

SHELTERS 

HUMBER  OF 
EXOATMOSFHERIC 
INTERCEPTORS 


Figure  S.  Layered  Experimental  Design 
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The  model  inputs  are  the  number  of  attacking  RVs,  the  number 
of  shelters,  and  the  number  of  exoatmospheric  interceptors. 
After  the  required  number  of  simulations,  the  model  outputs 
MX  survivability.  The  number  of  attacking  RVs  was  set  at 
either  16  (low  density  attack)  or  24  (high  density  attack). 
The  number  of  RVs  was  not  considered  a  design  factor  in  any 
of  the  experimental  designs  because  each  attack  level  was 
analyzed  separately.  The  number  of  shelters  was  set  at  a 
discrete  level  from  8  to  23 ,  and  the  number  of  interceptors 
was  set  at  any  integer  value  from  one  to  the  number  of 
attacking  RVs.  The  experimental  model  was  run  with  20  dif¬ 
ferent  combinations  of  shelters  and  interceptors,  and  the  MX 
survivability  of  these  combinations  was  recorded.  Regression 
analysis  was  then  used  to  develop  an  equation  which  fits 
these  20  data  points,  using  the  number  of  shelters  and  inter¬ 
ceptors  as  the  independent  variables,  and  MX  survivability  as 
the  dependent  variable.  The  general  form  of  the  regression 
equation  is: 


Y  •  bo  *  blXl  *  b2XZ  *  ’ 


•  *  b„*„ 
n  n 


where 

bQ  through  bR  ■  regression  coefficients 

Xj  through  Xn  *  a  form  of  the  independent  variables. 
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The  equation  which  characterizes  these  data  was  used  to 
determine  the  number  of  exoatmospheric  interceptors  needed 
to  obtain  a  particular  MX  survivability,  given  a  certain 
number  of  shelters  and  attacking  RVs.  A  form  of  this  equa¬ 
tion  was  selected  based  upon  its  large  adjusted  coefficient 
of  determination  (52),  since  the  main  concern  is  total 
predictive  power  and  not  the  marginal  predictive  power  of 
each  independent  variable.  A  regression  equation  with  an 
adjusted  R2  greater  than  approximately  0.64  is  considered 
to  have  significant  predictive  power  (Ref  39). 

Endoatmospheric  Defense  Model  (One  TDU) .  This  model’s 
experimental  design  consists  of  only  one  factor  (number  of 
shelters)  and  one  response  (see  Figure  6),  since  the  number 
of  TDUs  is  fixed  at  one.  The  inputs  to  the  model  are  the 
number  of  attacking  RVs  and  the  number  of  shelters  per  com¬ 
plex,  and  the  model  output  is  MX  survivability.  The  number 
of  RVs  is  set  at  either  16  (low  density  attack)  or  24  (high 
density  attack),  and  the  number  of  shelters  will  be  at  a 
discrete  level  from  8  to  23.  The  model  was  run  for  these 
different  shelter  values,  and  the  MX  survivability  for  these 
values  was  documented.  Regression  analysis  was  again  used 
to  develop  an  equation  which  fits  these  data  points,  with 
the  number  of  shelters  being  the  independent  variable  and  MX 

t 

i 

survivability  being  the  dependent  variable.  In  both  the 
one  and  two  TDU  cases,  where  all  possible  model  outcomes  are 
determined  by  simulation,  regression  analysis  provides  a  more 
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MX 

Survivability 


Number 
of  Shelters 


Figure  6.  One  TDU  Experimental  Model 

representative  output  by  smoothing  the  stochastic  processes 
of  the  model.  The  general  form  of  the  regression  equation 
is : 


Y  -  b0  .  Vl  ♦  b2X2  .....  bnXn  . 

A  form  of  this  equation  was  then  used  to  determine  the 
number  of  shelters,  which,  when  combined  with  the  one  TDU, 
provides  a  predetermined  level  of  MX  survivability. 

Endoatmospheric  Defense  Model  CTwo  TDUs).  The  experi¬ 
mental  design  of  the  two  TDU  model  consists  of  one  factor 
(number  of  shelters)  and  one  response  (MX  survivability)  as 
depicted  in  Figure  7.  The  model  inputs  are  the  number  of 
attacking  RVs,  which  is  set  at  16  (low  density  attack),  or 
24  (high  density  attack),  and  the  number  of  shelters  per 
complex,  which  is  set  at  a  discrete  value  from  8  to  23. 

The  model  output  is  MX  survivability.  As  in  the  case  of 
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Figure  7.  Two  TDU  Experimental  Design 


one  TDU,  the  model  was  executed  for  these  shelter  values, 
and  the  respective  MX  survivability  was  recorded.  Regression 
analysis  was  applied  to  these  data  points  to  develop  an 
equation  with  a  large  R2  ;  with  the  number  of  shelters 
being  the  independent  variable  and  MX  survivability  being 
the  dependent  variable.  Once  again,  the  general  form  of 
the  regression  equation  is: 


+  Vl  +  b2X2  + 


+  Mn 
n  n 


A  form  of  this  equation  was  employed  to  determine  the  number 
of  shelters,  which,  when  deployed  with  two  TDUs,  yields  a 
predetermined  level  of  MX  survivability  against  an  enemy 
attack. 
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Number  of  Sample  Runs 

The  effectiveness  of  all  the  models  is  dependent  upon 
the  number  of  runs  required  to  achieve  a  particular  level  of 
confidence  in  the  results.  The  results  of  the  models  can  be 
characterized  by  the  binomial  distribution  because  the  output 
of  one  model  run  is  a  Bernoulli  trial.  For  reasonably  large 
sample  sizes,  the  binomial  distribution  can  be  approximated 
by  the  normal  distribution.  It  can  be  shown  that,  in  a  worse 
case  situation,  the  equation  in  Chapter  II  (verification) 
reduces  to: 

N  *  Z^/2/C4e2) 

where  N  is  the  maximum  number  of  model  runs  and  is  an  upper 
bound  for  any  degree  of  confidence,  e  is  the  desired  inter¬ 
val  (in  decimal)  about  the  true  mean,  and  ^a/2  is 
standardized  normal  statistic  for  the  probability  sought 
(Ref  32:191-192).  Since  it  is  desired  that  the  model  results 
differ  from  the  true  mean  by  less  than  approximately  3%  with 
a  95%  confidence  interval  (za/2  *  1-96)  ,  the  maximum  number 
of  model  runs  required  is  1067  or  approximately  1200. 

Cost  Effectiveness 

The  following  cost  data  (Ref  1:20)  will  be  applied  to 
the  parameters  of  the  layered  and  endoatmospheric  defense 
systems  which,  according  to  the  regression  equations  selected, 
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provide  a  predetermined  level  of  MX  survivability. 


Applying  this  cost  data  results  in  a  total  cost  figure 
for  each  defense  configuration.  A  comparison  of  these  costs 
will  provide  the  most  cost-effective  solution  for  defending 
the  MX,  since  all  of  the  alternatives  considered  provide  an 
equivalent  level  of  survivability. 


TABLE  VII 

Costs  of  Offensive  and  Defensive  Systems 


Parameter 

Fixed  Cost 
($M) 

Variable  Cost 
($M) 

MX  Missile 

8000 

60  (M)  ‘  7  8 

Shelter 

5000 

3S 

Exoatmospheric 

Interceptor 

7000 

.78 

60  (X) 

Endoatmospheric 

Interceptor 

5000 

16(N)*78 

where 

M  *  number  of  MX  missiles 
S  *  number  of  shelters 

X  -  number  of  exoatmospheric  interceptors 
N  ■  number  of  endoatmospheric  interceptors  (Ref  1:20). 
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Model  Runs 


The  layered  defense  model  was  run  1200  times  for  twenty 
different  combinations  of  shelters  and  exoatmospheric  inter¬ 
ceptors  being  attacked  by  16  and  24  Soviet  RVs.  For  both  the 
one  and  two  TDU  endoatmospheric  models,  the  model  was  exe¬ 
cuted  1200  times  for  the  twelve  shelter  configurations  (12 
through  23)  under  attack  by  24  Soviet  RVs,  and  1200  times  for 
the  16  shelter  configurations  (8  through  23)  under  attack  by 
16  Soviet  RVs.  The  lower  shelter  value  is  restricted  due  to 
the  fratracide  limit  of  a  maximum  of  two  RVs  being  targeted 
against  each  shelter. 

Results 

The  output  of  a  model  simulation  is  MX  survivability 
(MXS) ,  and  the  results  of  the  model  runs  mentioned  above 
are  presented  in  Tables  VIII  through  XIII. 


[Layered  Defense) 


Number  of 
Shelters 

Number  of 

Exoatmospheric  Interceptors 

MX  Survivability 
(%) 

12 

1 

44.3 

12 

4 

48.9 

12 

5 

51.8 

12 

6 

56.1 

12 

8 

61.5 

13 

1 

57.9 

13 

2 

62.7 

13 

15 

82.8 

14 

1 

61.4 

14 

10 

73.7 

15 

11 

75.3 

16 

20 

83.6 

17 

4 

68.0 

17 

16 

80.8 

17 

23 

87.7 

18 

20 

84.5 

20 

20 

84.5 

23 

20 

85.4 

23 

23 

94.0 

23 

24 

99.1 
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TABLE  IX 

MX  Survivability  Against  a  24  RV  Attack 
(One  TDU  Defense) 


Number  of 

Shelters 

MX  Survivability 
(%) 

12 

40.9 

13 

57.4 

14 

58.8 

15 

60.1 

16 

60.3 

17 

60.0 

18 

63.0 

19 

62.6 

20 

62.1 

21 

63.0 

22 

64.2 

23 

71.9 

9 


TABLE  X 

MX  Survivability  Against  a  24  RV  Attack 
(Two  TDU  Defense) 


Number  of 

Shelters 

MX  Survivability 

m 

12 

50.3 

13 

60.0 

14 

65.3 

15 

65.0 

16 

68.0 

17 

65.0 

18 

66.7 

19 

71.3 

20 

69.2 

21 

68.8 

22 

69.7 

23 

73.3 
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TABLE  XI 

MX 

Survivability  Against  a  16  RV 

Attack 

(Layered  Defense) 

Number  of 
Shelters 

Number  of 

Exoatmospheric  Interceptors 

MX  Survivability 
(*) 

8 

1 

45.4 

8 

2 

49.8 

8 

5 

59.7 

9 

1 

62.5 

10 

2 

63.5 

10 

10 

80.5 

12 

S 

70.4 

12 

16 

92.5 

14 

1 

63.4 

IS 

12 

85.8 

16 

8 

83.5 

17 

4 

87.3 

17 

12 

93.7 

18 

10 

93.2 

20 

5 

92.0 

20 

14 

96.8 

20 

16 

98.0 

23 

2 

89.4 

23 

8 

94.2 

23 


16 


97.4 


TABLE  XII 

MX  Survivability  Against  a  16  RV  Attack 
(One  TDU  Defense) 


Number  of 

Shelters 

MX  Survivability 
(*) 

8 

40.4 

9 

60.7 

10 

58.1 

11 

60.9 

12 

62.2 

13 

62.0 

14 

64.5 

15 

69.3 

16 

68.9 

17 

69.9 

18 

85.1 

19 

84.0 

20 

84.7 

21 

85.6 

22 

86.6 

23 

86.7 
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TABLE  XIII 

MX  Survivability  Against  a  16  RV  Attack 
(Two  TDU  Defense) 


Number  of 

Shelters 

MX  Survivability 
(%) 

8 

40.4 

9 

62.6 

10 

63.7 

11 

67.0 

12 

67.4 

13 

70.4 

14 

68.4 

15 

74.0 

16 

76.5 

17 

80.9 

18 

84.9 

19 

86.3 

20 

85.6 

21 

87.1 

22 

87.6 

23 

87.7 

Analysis 

The  experimental  models  discussed  previously  were  used 
to  develop  the  equations  which  characterize  the  results  for 
each  of  the  defense  systems.  Various  equations  Csemilog, 
reciprocal,  polynomial,  and  linear)  which  might  characterize 
the  true  relationship  of  the  data  were  evaluated.  Regression 
analysis  was  used  because  the  equations  obtained  from  this 
process  provide  the  following  advantages: 

1.  Conserves  computer  resources  (i.e.,  limits  the 
number  of  model  runs  required). 

2.  Smooths  the  stochastic  processes  of  the  models. 

3.  Easily  applicable  to  different  survivability 
levels  without  further  computer  simulations. 

The  layered  defense  equation  can  be  used  to  develop  the 
combinations  of  the  number  of  exoatmospheric  interceptors  and 
shelters  which  provide  a  prescribed  level  of  survivability, 
while  both  the  one  and  two  TDU  equations  can  be  used  directly 
to  obtain  the  number  of  shelters  which  provide  the  prescribed 
survivability.  For  assured  destruction,  the  U.S.  retaliatory 
strike  force  has  been  arbitrarily  chosen  to  consist  of  approx¬ 
imately  1200  warheads  that  can  be  delivered  against  Soviet 
targets  of  value,  or  120  MX-equivalent  payloads  (Ref  1:10). 
This  equates  to  a  prescribed  MX  survivability  of  601.  When 
using  the  equations  to  determine  the  number  of  interceptors 
and/or  shelters  which  provide  60%  MX  survivability,  the 
result  may  consist  of  a  fraction  of  an  interceptor  or  shelter. 
Since  60%  is  the  chosen  requirement  for  U.S.  MX  survivability, 
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the  general  rule  for  rounding  off  in  this  study  is  to  round 
up  to  the  next  highest  number  of  interceptors  or  shelters. 
Fractions  of  interceptors  or  shelters  less  than  one-tenth 
are  assumed  to  be  negligible  and  therefore  will  be  rounded 
down  to  the  nearest  integer  value. 

High  Density  Attack  (24  RVs). 

Layered  Defense  Model.  A  regression  analysis  of 
the  20  data  points  collected  for  this  model  was  performed 
using  the  Statistical  Package  for  the  Social  Sciences  (SPSS). 
Of  the  equations  evaluated,  the  following  model  best  charac¬ 
terizes  the  true  relationship  of  the  data: 


MX  Survivability  -  bQ  +  +  b2X2  +  **3*3 *  R2  *  0.905 


where 


b 

b 

b 

b 

X 

X 

X 


o 

1 

2 

3 

1 

2 

3 


1.034 

.02173 

-6.518 

-.000632 


F-Ratio 

44.57 

13.46 

10.90 

3.14 


number  of  exoatmospheric  interceptors 
1/number  of  shelters 
X^*  (number  of  shelters) 


This  equation  was  selected  because  it  has  the  best  predictive 
power  (i.e.,  largest  significant  adjusted  R2)  and  all  of  the 
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parameters  estimated  by  this  model  (bQ,  b^,  b2»  and  bj)  were 
significantly  different  than  zero  (at  90%  confidence)  as 
indicated  by  their  partial  F-ratios. 

Using  the  equation  selected.  Table  XIV  shows  the  combin¬ 
ations  of  exoatmospheric  interceptors  and  shelters  which 
provide  at  least  60%  MX  survivability.  All  of  the  combina¬ 
tions  which  do  not  require  any  exoatmospheric  interceptors 
to  achieve  60%  MX  survivability  are  actually  an  endoatmos- 
pheric  defense  with  one  TDU  and  not  a  layered  defense. 


TABEL  XIV 

60%  MX  Survivability 
(Layered  Defense  -  24  RV  Attack) 

Number  of 

Number  of 

Shelters 

Exoatmospheric  Interceptors 

12 

8 

13 

5 

14 

3 

15 

0 

Therefore,  for  the  layered  defense  model,  the  cost  data  will 
only  be  applied  to  the  three  combinations  which  require  exo¬ 
atmospheric  interceptors.  Table  XV  presents  the  total  cost 
to  provide  all  200  MX  complexes  with  any  of  the  three  layered 
defenses. 
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TABLE  XV 

Total  Cost  of  200  MX  Complexes 
(Layered  Defense  -  24  RV  Attack) 


Number  of 
Shelters 

Number  of 

Exoatmospheric  Interceptors 

Total  Cost 
($  Millions) 

12 

8 

57,230 

13 

5 

52,017 

14 

3 

43,303 

This  table  shows  that  a  defense  consisting  of  14  shelters, 
three  exoatmospheric  interceptors,  one  TDU,  and  one  MX  in 
each  of  the  200  MX  complexes  is  the  most  cost-effective 
layered  defense  option  against  a  24  RV  attack. 

Endoatmospheric  Defense  Model  (One  TDU).  An 
equation  which  fits  the  twelve  data  points  was  obtained 
using  SPSS,  and  the  following  equation  was  chosen  as  the 
best  fit  for  the  data: 

MX  Survivability  *  b^  +  +  ^2*2  +  ^3*3  >  ■  0.823 

where  F-Ratio 


bO  * 

11.0566 

26.6 

bl  * 

-67.9109 

26.4 

b2  * 

-.44839 

22.1 

N  - 

.00024 

20.7 
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1/number  of  shelters 


number  of  shelters 
number  of  shelters  cubed 

Of  all  the  models  examined,  the  equation  selected  has  the 
bestT  predictive  power  (i.e.,  largest  significant  K2)  and  the 
partial  F-ratios  of  this  model  indicate  that  all  of  the  model 
parameters  estimated  (i.e.,  bQ,  b^,  b2,  and  b3)  were  signi¬ 
ficantly  different  from  zero. 

Substituting  60%  MX  survivability  into  this  equation 
yields  14.39  shelters.  Therefore,  15  shelters  are  needed 
to  obtain  60%  MX  survivability  with  one  terminal  defense 
unit.  The  total  cost  for  developing  all  200  MX  complexes 


with  one  TDU 

presented  in 

per  complex  and  15  shelters  per  complex  is 

Table  XVI. 

TABLE  XVI 

Total  Cost  of  200  MX  Complexes 

(One  TDU  Defense  -  24  RV  Attack) 

Number  of 
Shelters 

Number  of 

Endoatmospheric  Interceptors 

Total  Cost 
($  Millions) 

15 

3 

33,090.77 
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Endoatmospheric  Defense  Model  (Two  TDlTs).  Applying 
regression  analysis  to  the  twelve  data  points  using  SPSS 
resulted  in  the  following  equation  having  the  best  adjusted 
coefficient  of  determination  (R2)  and  model  parameters  (i.e., 
b  ,  b^,  b2»  and  bj)  significantly  different  from  zero: 

MX  Survivability  *  bQ  +  b^X^  +  ^2^2  +  b3^3  *  *  0-897 


where 

F -Ratio 

bo 

a 

6.5358 

12.9 

bl 

a 

-40.0277 

12.68 

b2 

a 

-.24153 

8.86 

b3 

a 

.0001225 

7.45 

X1 

a 

1/number  i 

of  shelters 

X2 

a 

number  of 

shelters 

X3 

a 

number  of 

shelters  cubed. 

Other  models  examined  were  semilog,  linear,  and  log  models. 

For  60%  MX  survivability,  the  equation  selected  yields 
13.26  shelters.  Hence,  14  shelters  are  required  to  obtain 
60%  MX  survivability  with  two  terminal  defense  units. 

Table  XVII  shows  the  total  cost  for  defending  all  200  MX 
complexes  with  two  TDUs  per  complex  and  14  shelters  per 
complex . 


TABLE  XVII 

Total  Cost  of  200  MX  Complexes 
(Two  TDU  Defense  -  24  RV  Attack) 


Number  of 

Number  of 

Total  Cost 

Shelters 

Endoatmospheric  Interceptors 

($  Millions) 

14 

6 

34,176.07 

Low  Density  Attack  (16  RVs). 

Layered  Defense  Model.  A  regression  analysis  of 
the  20  data  points  collected  for  this  model  was  performed 
using  SPSS.  Of  the  equations  evaluated,  the  following  model 
best  characterizes  the  true  relationship  of  the  data: 

MX  Survivability  -  bQ  +  b.^  +  b2X2  +  b^Xj  ,  TT2  »  0.937 


where 

F- Ratio 

bo 

m 

.  24729 

31.5 

bl 

m 

.02974 

92.78 

b2 

M 

.03971 

40.03 

b3 

m 

-.0015589 

17.57 

X1 

m 

number  of 

shelters 

X2 

m 

number  of 

exoatmospheric  interceptors 

X3 

m 

60 


This  equation  was  selected  because  it  had  the  best  predictive 
power  and  all  of  the  parameters  estimated  by  this  model  (i.e., 
b0>  b^,  b2>  and  bj)  were  significantly  different  than  zero 
(at  90%)  as  indicated  by  their  partial  F-ratios. 

Using  the  equation  selected,  Table  XVIII  shows  the 
combinations  of  exoatmospheric  interceptors  and  shelters 
which  provide  60%  MX  survivability.  All  of  the  combinations 
which  do  not  require  any  exoatmospheric  interceptors  to 
achieve  60%  MX  survivability  are  again  actually  an  endo- 
atmospheric  defense  with  one  TDU  and  not  a  layered  defense. 


TABLE  XVIII 

60%  MX  Survivability 
(Layered  Defense  -  16  RV  Attack) 


Number  of  Number  of 

Shelters  Exoatmospheric  Interceptors 


8  5 

9  4 

10  3 

11  2 

12  0 


Therefore,  for  the  layered  defense  model,  the  cost  data  will 
only  be  applied  to  the  four  combinations  which  require  exo¬ 
atmospheric  interceptors.  Table  XIX  presents  the  total  cost 
to  provide  all  200  MX  complexes  with  any  of  the  layered  defenses. 
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TABLE  XIX 


Total  Cost  of  200  MX  Complexes 


Number  of 
Shelters 

Number  of 

Exoatmospheric  Interceptors 

Total  Cost 
($  Millions) 

8 

5 

49,017 

9 

4 

47,520 

10 

3 

45,903 

11 

2 

44,114 

This  table  shows  that  a  defense  consisting  of  11  shelters 
two  exoatmospheric  interceptors,  one  TDU,  and  one  MX  in 
each  of  the  200  MX  complexes  is  the  most  cost-effective 
layered  defense  option  against  a  16  RV  attack. 

Endoatmosoheric  Defense  Model  fOne  TDU') .  An 
equation  which  fits  the  16  data  points  was  obtained  using 
SPSS,  and  the  following  equation  was  chosen  as  the  best 
fit  for  the  data: 

MX  Survivability  ■  bQ  +  b^X^  ,  IT2  *  0.881 


where 


F- Ratio 


b  -  .2884  49.26 

o 


.02695  112.46 


X, 


number  of  shelters 


Of  all  the  models  examined,  the  equation  selected  has  the 
best  predictive  power  and  the  partial  F-ratios  of  this  model 


indicate  that  all  of  the  model  parameters  estimated  (i.e., 
bQ  and  b^)  are  significantly  different  from  zero. 

Substituting  60%  MX  survivability  into  this  equation 
yields  11.56  shelters.  Therefore,  12  shelters  are  needed 
to  obtain  60%  MX  survivability  with  one  terminal  defense 
unit.  The  total  cost  for  defending  all  200  MX  complexes 
with  one  TDU  per  complex  and  12  shelters  per  complex  is 
presented  in  Table  XX. 


TABLE  XX 

Total  Cost  of  200  MX  Complexes 
(One  TDU  Defense  -  16  RV  Attack] 


Number  of 

Number  of 

Total  Cost 

Shelters 

Endoatmospheric  Interceptors 

($  Millions) 

12 

3 

31,290.77 

Endoatmospheric  Defense  Model  (Two  TDUs).  Applying 
regression  analysis  to  the  16  data  points  using  SPSS  resulted 
in  the  following  equation  having  the  best  adjusted  coefficient 
of  determination  (E2)  and  model  parameters  (i.e.,  bQ  and  b^) 
significantly  different  from  zero  at  the  901  level: 
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MX  Survivability  ■  bQ  +  bjX^^ 


9 


R2  -  0.948 


where 


F-Ratio 


-.151528 


7.613 


,334915 


273.42 


*  An  (shelters). 


For  60$  MX  survivability,  the  equation  selected  yields 
9.43  shelters.  Hence,  10  shelters  are  required  to  obtain 
60$  MX  survivability  with  two  terminal  defense  units. 

Table  XXI  shows  the  total  cost  for  defending  all  200  MX 
complexes  with  two  TDUs  per  complex  and  10  shelters  per 
complex. 


TABLE  XXI 


Total  Cost  of  200  MX  Complexes 
(Two  TDU  Defense  -  16  RV  Attack) 


Number  of 

Number  of 

Total  Cost 

Shelters 

Endoatmospheric  Interceptors 

($  Millions) 

10 

6 

31,776.07 

6 


Summary 

High  Density  Attack  (24  RVs).  Of  the  three  defense 
systems  investigated,  the  one  TDU  defense  system  is  the 
most  cost-effective.  The  terminal  BMD  system  with  two  TDUs 
is  almost  as  cost-effective  as  one  TDU,  but  both  terminal 
BMD  strategies  are  much  more  cost  effective  than  the  layered 
defense  strategy. 

Low  Density  Attack  (16  RVs).  Of  the  three  defense 
systems  investigated,  the  one  TDU  defense  system  is  again 
the  most  cost-effective  and,  as  in  the  high  density  attack, 
both  terminal  BMD  strategies  are  much  more  cost-effective 
than  the  layered  defense  strategy. 


IV.  Conclusions 


Based  on  the  results  obtained  from  this  study,  the 
development  cost  of  the  exoatmospheric  component  causes  the 
layered  defense  system  to  be  non-competitive  with  the  term¬ 
inal  defense  systems  when  cost  is  the  only  criteria.  It 
is  possible  that  a  substantial  portion  of  the  development 
cost  could  be  overcome  if  the  exoatmospheric  defense  com¬ 
ponent  could  preferentially  defend  the  MX  complex  (i.e., 
intercept  only  those  RVs  aimed  at  itself  or  the  actual  MX) . 
If  preferential  defense  by  the  exoatmospheric  defense  layer 
can  be  achieved,  the  end oatmo spheric  portion  of  the  layered 
defense  can  be  eliminated  because  of  the  relatively  high  PK 
of  an  exoatmospheric  interceptor  with  a  three  foot  effective 
weapons  radius  and  one  and  one-half  foot  CEP.  Hence,  the 
layered  defense  system  would  become  solely  an  exoatmospheric 
defense  system. 

The  two  TDU  effectiveness  (i.e.,  MX  survivability)  was 
only  marginally  better  than  the  one  TDU  effectiveness,  pri¬ 
marily  because  there  are  three  shelters  in  each  complex  to 
defend  rather  than  two  as  in  the  one  TDU  case.  Additional 
factors  to  consider  here  are  the  possibility  of  detonating 
twice  as  many  nuclear  weapons  (six  vice  three)  at  a  low 
altitude  over  friendly  territory  with  two  TDUs,  and  the 
increased  amount  of  high-neutron  yielding  fissile  material 
that  must  be  withdrawn  from  offensive  warhead  production. 


Due  to  the  nuclear  kill  mechanism  of  the  LoADs  system, 
it  is  only  effective  in  defending  hard  targets  and  cannot 
be  used  to  defend  soft  area  targets.  The  exoatmospheric 
interceptors  help  in  this  regard  and  can  reduce  the  number 
of  nuclear  atmospheric  detonations  required.  A  superior 
system  would  preferentially  intercept  all  RVs  at  high  altitudes 
without  nuclear  kill,  thus  being  able  to  defend  soft  area 
targets  as  well  as  hard  point  targets.  The  results  do  clearly 
indicate  that  significant  flexibility  exists  for  tradeoff 
between  the  number  of  protective  shelters  and  various  defense 
configurations  for  a  preset  level  of  MX  survivability. 

While  the  actual  number  of  shelters  or  interceptors 
generated  by  the  regression  models  in  this  study  provide 
an  indication  of  the  magnitudes  involved,  these  regression 
model  outputs  vary  from  the  simulation  model  outputs  and, 
as  with  all  model  outputs,  should  be  employed  with  caution. 
Therefore,  absolute  recommendations  should  not  be  made,  but 
the  relative  comparisons  of  the  layered,  one  TDU,  and  two 
TDU  defense  systems  in  this  study  can  be  drawn. 

Several  very  sensitive  variables  were  discovered  in  the 
course  of  this  study  and  the  values  used  for  these  variables 
will  determine  the  outcome  of  the  models.  These  variables 
were  sure-safe  and  sure-kill  neutron  intensity  levels  of  the 
RVs,  CEPs  for  the  RVs  and  all  of  the  interceptors,  and  the 
effective  weapon  radius  of  the  exoatmospheric  interceptors. 
Shelter  hardness  was  not  a  critical  variable  as  long  as  the 
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sure-safe  and  sure-kill  intensities  remained  500  psi  apart. 
Values  closer  together  would  make  this  variable  more 
sensitive. 

While  the  layered  defense  system  is  not  justified  on 
cost  considerations  alone,  the  research  required  in  building 
the  system  could  lead  to  tremendous  improvements  such  as 
preferential  defense  by  the  exoatmospheric  defense  layer, 
direct  impact /conventional  warhead  kill  within  the  atmosphere, 
and  an  absolute  degree  of  confidence  in  the  system's  ability 
because  it  could  be  tested  in  actual  RV  intercept  situations 
where  the  nuclear  intercept  cannot. 
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V.  Recommendations  for  Further  Study 


Soviet  Union  parity  with  the  U.S  in  the  nuclear  arms  race 
has  stimulated  considerable  interest  in  BMD  as  an  addition  to 
the  current  MPS  basing  mode;  and  hence  efforts  in  BMD,  beyond 
this  study,  could  prove  very  helpful.  A  scenario  using  clas¬ 
sified  values  for  the  parameters  and  variables  of  these  models 
is  a  necessity  for  an  accurate  representation  of  reality. 

Various  areas  of  the  analysis  and  models  in  this  study 
which  could  be  expanded  or  enhanced  are  stated  below: 

1.  The  RV-shelter  interaction  of  the  model  could  be 
enhanced  to  include  ground  shock  as  another  possible 
RV-kill  mechanism. 

2.  The  survivability  and  attack  levels  chosen  in  the 
analysis  phase  could  be  expanded  to  investigate 
several  attack  and  defense  scenarios. 

3.  The  cost  data  used  in  the  analysis  phase  could  be 
enhanced  by  obtaining  more  detailed  cost  estimates. 

4.  The  computer  model  could  be  enhanced  to  more 
efficiently  assign  RVs  to  targets  for  all  attack 
and  defense  configurations.  This  would  eliminate 
the  model  modifications  in  Appendix  G. 

5.  The  computer  model  could  be  expanded  to  include 
other  possible  interceptor  strategies. 

Finally,  as  both  the  exoatmospheric  and  endoatmospheric 
defense  subsystems  are  researched,  developed,  and  refined, 
the  requirements  for  future  efforts  in  the  area  of  BMD  will 
expand . 
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Appendix  A 


Probability  of  Kill  Due  to  Cratering 


MX  Shelter 

A  surface  burst  nuclear  explosion  creates  a  crater. 

The  crater  ejecta  consists  of  soil  or  rock  debris  that  are 
thrown  beyond  the  boundaries  of  the  apparent  crater  (Ref  18: 
255).  It  is  assumed  that  the  shelter  door  shown  in  Figure  8 
will  be  totally  covered  at  2/3  the  radius  of  the  ejecta  (Re) 
since  ejecta  thickness  is  approximately  16.86  feet  at  0.84  Re. 
In  addition,  any  object  within  the  apparent  crater  radius 
(R  )  will  be  destroyed. 

To  calculate  both  R„  and  R  ,  the  following  equations 

6  cl 

are  used: 


Ra  -  RaRCr3) 

Re  -  2.15  Ra 

where  Y  is  the  yield  of  the  weapon  in  kilotons  and  Rfl 

is  the  apparent  crater  radius  created  by  a  one  kiloton  weapon. 

R 

For  a  surface  burst  in  dry  soil  or  dry  hard  rock,  Ra  is 
61  feet  (Ref  18:254-256).  The  probability  of  kill  (PK)  for 
cratering  is  defined  by  the  circular  normal  function  and  can 
be  found  using  the  following  equation. 


Mia3aK.g  ft 3 
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+  (Ra+85.S) 


PK 


/ 


M-Z-y 


+  (R +66.5) 

2L 


JZn  a 

(2/3Re+85.S)  y 


dy 


/ 


/TfF  a 


e  x  dx 


-(R  +66.  S) 

cl 


where 

a  *  o  *  CEP  /2  Jin  2  *  a 

x  y 

Letting  zi  *  ^  ■+  dx  *  adz  ,  and  z2  *  £  ■+  dy  ■  adz 
and  substituting  into  the  above  equation  yields 


+  (R  +85.5)  +CR+66.5) 

8  cL 


-(2/3R+85.5)  -(R  +66.5) 

_  c  21 


a  a 


The  MX  integration  limits  in  the  y  direction  were  chosen 

assuming  the  apparent  crater  radius  (R_)  had  to  at  least 

touch  the  back  of  the  shelter  for  any  chance  of  destruction 

from  cratering  (+y  direction)  and  the  ejecta  at  2/3  Rg 

has  to  cover  the  door  (-y  direction).  For  the  x-direction, 

it  was  assumed  that  the  R„  must  at  least  touch  the  sides 

a 

of  the  shelter. 

Using  a  circular  error  probable  (CEP)  of  0.2  NM 
(1215.2  feet),  a  yield  of  1000  kilotons,  and  the  dimensions 
of  an  MX  shelter  (Figure  8) ,  the  following  results  are 
obtained: 


5 


R  -  (61  feet)  (1000  )  -  484.54  feet 

cl 

Re  *  (494.54  feet)  (2.15)  -  1041.76  feet 

2/3Re  »  694.51  feet 

o_  -  Q.  *  121.5:2  flet.  *  1032.1  feet.  • 

x  V  SI  "in"  2 


Using  a  standard  normal  table  (Ref  37), 

PK  -  (.7096  -  .2248)  X  (.7033  -  .2967) 


PK  -  (.4848)  X  (.4066) 


.1971. 


Terminal  Defense  Unit  (TDU)  Shelter 


A  TDU  shelter  is  susceptible  to  the  same  cratering 
effects  as  the  MX  shelter,  except  that  the  limits  of  inte¬ 
gration  for  cratering  a  TDU  shelter  will  be  different 
because  the  TDU  is  designed  to  punch  out  of  the  top  of  an 
MX  shelter  (Figure  9)  (Ref  37).  Covering  any  portion  of 


Figure  9.  TDU  Shelter  (Ref  37) 

the  shelter  with  ejecta  is,  therefore,  assumed  to  be  inef¬ 
fective  in  prohibiting  TDU  operation.  The  probability  of 


78 


— I|TWW> 


kill  (PK)  of  cratering  is  defined  by  the  circular  normal 
function  and  can  be  found  using  the  following  equation: 


+  (R  +85.5)  +(R  +66,5) 

a  a 


-(R+8S.5)  -(R+66.5) 

a  a 


a  a 

The  TDU  integration  limits,  in  both  directions  were  chosen 
assuming  R  must  touch  the  perimeter  of  the  shelter  for 

a 

any  chance  of  destruction  from  cratering. 

Using  the  same  CEP,  yield,  R  ,  and  R  as  the  MX 
calculations,  the  probability  of  killing  (PK)  a  TDU  can  be 
calculated  as  follows: 
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PK 


+.5523 


-.5523 


+  .534 


Using  a  standard  normal  table  (Ref  37), 

PK  -  (.7096  -  .2904)  X  (.7033  -  .2967) 

PK  -  (.4192)  X  (.4066)  -  .1704 
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Appendix  B 


Ten  Cell  Model 

A  methodology  used  to  calculate  the  probability  of 
kill  (PK)  for  a  particular  weapon  against  a  specific  target 
is  based  on  the  ten  cell  model.  This  model  places  ten  cells 
of  equal  probability  of  hit  centered  around  the  designated 
ground  zero  (DGZ) ,  which  is  the  aim  point  of  the  weapon. 
Thus,  the  weapon  has  a  10%  chance  of  impacting  each  cell. 

The  inputs  required  for  this  model  are  the  sure- safe  and 
sure-kill  weapon  effect  intensities  of  the  target,  the  CEP 
of  the  weapon,  and  the  distance  of  the  target  of  interest 
from  the  DGZ.  Figure  10  depicts  this  cell  structure. 


Figure  10.  Ten  Cells  of  Equal  Probability  of  Hit 


Note  that  the  distance  from  DGZ  is  infinity  for  the  outer¬ 
most  circle  of  the  model.  The  case  where  a  target  is  a 
given  distance  from  the  DGZ  and  a  weapon  impacts  at  some 
third  point  is  the  most  general  case  and  is  depicted  in 
Figure  11.  The  PK  of  a  target  at  a  distance  X  from  the 


DGZ  is  a  function  of  the  probability  of  hit  (PH)  at  a 
point  described  by  p  ,  <p  ,  and  the  probability  of  damage 
(PD)  of  the  target  at  a  distance  r  from  the  point  of 
impact.  Thus,  the  following  equation  describes  the  situation: 

PK(X)  -  (Prob  of  Hitting  Cell  i) 

X  (Prob  of  Damage  Given  a  Hit  on  Cell  i) 


or 
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2ir  * 

PK(X)  -  /  /  PH(pt$)  PDCr)rdrd<|> 

0  0 


Placing  the  ten  cell  model  over  Figure  11  and  treating  each 
cell  as  a  discrete  impact  point  produces  the  following 
equation: 


PK(X) 


T 

l  PH(P. AA-  PDfr.) 
i*l  iix  i 


CD 


If  we  know  a  probability  damage  function  based  on  intensity 
(PD (I))  ,  this  can  be  substituted  directly  for  the  proba¬ 
bility  damage  function  based  on  range  (PD(r))  .  The 
variables  are  defined  as  follows: 

r^  ■  distance  from  target  to  the  center  of  cell  i; 

I  ■  weapon  effect  intensity  at  the  center  of  cell  i; 

p-  -  distance  from  the  DGZ  to  the  outer  edge  of 
cell  i; 

<J>^  *  angle  of  cell  i  in  relation  to  the  DGZ; 

AA^  *  area  of  cell  i; 

Nt  »  number  of  cells  in  the  model; 

<p.>  *  distance  from  the  DGZ  to  the  probabilistic 

center  of  cell  i; 

<$ .>  *  angle  at  which  the  probabilistic  center  of 

cell  i  is  located; 

n^  *  number  of  cells  in  ring  i; 

N.  •  number  of  cells  in  ring  i  plus  all  cells 

A  inside  of  ring  i. 

Figure  12  illustrates  the  geometry  of  these  variables. 
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Figure  12.  Ten  Cell  Model  on  a  Target 
For  the  ten  cell  model,  the  following  values  can  be  assigned 

NT  -  10  Nx  -  1  N2  -  S  N3  -  10 

-  1  n2  »  4  n3  »  5  P3  ■  00 
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PH(p^,<|k) ■  1/N^.  for  each  i  since  the  model  is 
constructed  so  that  an  attacking  weapon  has  an  equal 
probability  of  hitting  within  each  cell.  Thus,  Eq  (1) 
can  be  written: 

1  Nt 

PK(X)  “  I  PD(I)  .  From  the  Law  of  Cosines, 

T  i-1 

r?  *  <p^>2  +  x2  -  2  <pi>  x  cos  <<p^> 


To  solve  for  <p^>  ,  pi  must  be  found.  Since  each  cell 
of  the  model  represents  an  equal  probability  of  hit  and 
since  the  probability  of  hit  is  distributed  circular  normal, 
the  following  equality  can  be  used: 


P(hit  from  DGZ  to  Ring  i) 


N. 

l 

Rip 


2  7T  pi 

/  / 

0  0 


.-Mp /o)! 
2tto2 


pdpde 


Integrating  over  9  and  simplifying  gives: 


N. 

x 

FTjT 


-l  / 1  t-Moy  pdP 

a2  0 


Letting  Z  ■  p/a  (implies  dZ  ■  dp/a)  and  solving  for 
PA  gives: 
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N. 


l 


N-p 


P  Jo 

-  /  e'5*2 2  (-Z )dZ 

0 


-^(p-./a)2 
-Ce  1 

-^(Pi/a) 

1  -  e  1 


e 


0 

3 


2 


Simplifying  results  in  the  following  expression  for  pi  : 

•h(p./c)2  N  N. 

e  or  (P^a)2  *  -2*n  (1  -  .  (2) 

To  find  an  expression  for  o2  ,  we  note  that  the  CEP  can 
be  defined  as: 

4  >  ;EP  e-Mr/o)2  _r  dT  .  1  .  e->stCEP/o)‘ 

0  cr2 

Simplifying: 

(CEP/a)2  =  -2  (£nl  -  Jln2)  or  (CEP/a)2  -  2  Jln2  , 


which  implies 


a 


2 


CEP2 

TXnT 


(3) 


Substituting  this  last  result  into  Eq  (2)  and  simplifying 
yields : 
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AO-All*  694 
UNCLASSIFIED 


AIR  FORCE  INST  OF  TECH  WRIBHT-PATTERSON  AFB  OH  SCHOO— ETC  F/S  15/3 
MX  SURVIVABILITY*  PASSIVE  AND  ACTIVE  DEFENSE. IU) 

MAR  BE  E  F  RETTAMMEL*  J  F  SHEEDY 

AFIT/BST/OS/B1M-11  NL 


(4) 


-AnCl  -  N./Nt)  h 

C - — -) 

An  2 


Now  <p^  must  be  found  in  terms  of  known  quantities.  We 
are  finding  the  center  of  a  two-dimensional  Gaussian  curve  so: 


<Pi> 


/ 

pi-l 

pi 

pi-l 


p  TOT  e-^‘2,pdp 


^e'^P/OS.pdp 


the  denominator  is  just  the  fraction  of  cells  between 
rings  i-1  and  i  :  ni/N<r  .  Therefore, 


N_  pi 


<pt>  -  ?  I  '■%Cb/3’!  <-p  • 

i  Pi_x  «2 


Letting  Z  *  p/a  (therefore  dZ  *  dp /a)  and  when 
P  *  P±  f  Z  M  P ±/o 

nt  ^i/o  -1,1* 

<pi>  •  I  1  e  d*  • 

1  pi-l 


-kz2 

Integrating  by  parts  with  U  -  Z  and  dV  ■  e  ^  zdz 

-kZ21 

(implies  dU  ■  dZ  and  V  ■  -e  1  gives: 


<Pi> 


N  ,  pi/a 

£  a(C-Ze-«!3 

1  pi-l 


Pi/P 

/ 

Pi-1 


-  e'hZ2  dZ} 
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The  integral  in  Eq  (5)  can  be  written  as  the  difference 
between  cumulative  normal  functions: 


Pi/a 


’i-l 


-  /2tt 


Pi/o 

/ 

0 


/Tk 


dZ 


pi-l/a  -hi2 
J  ® -  dZ 

0  /lirr 


Therefore, 


<p±> 


Eq  (5)  can  be  rewritten  as 


JL  e  1  1 

n^  o 


pi  CPi/a3  2 
~  e 


♦  /E?  CF (—)  - 

where  F  represents  the  cumulative  normal  function.  From 

Eq  (3), 


CEP 

/Zin2 


and  substituting  into  the  last  equation,  yields 
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N„ 


<p  >  .  ll  _CM_  {rttt 

1  ni 


n"  /2Tn7 


CEF"  pi-l  e 


-£n2(pi.1/CEP) 


Simplifying: 


/THT  -2.n2(p./CEP): 
TIT"  pi  e 


_  *'2£n2  p  .  j-m  i 

+  /lilF  ( - -  FC/n51 


CEF"  pi-l 


P,-_  i)3> 


<p.> 

"Cep" 


NT  t Pi~l 

HT  *  CEP  e 


-£n2 (pi_1/CEP) : 


Pi  -&n2 (pj/CEP) 2 

CEF  e 


+  /tt/ in 2  U  (^|pZ|  P^)  '  F  pi-l^^  (6) 


Because  the  two-dimensional  Gaussian  distribution  is  symmetric 
with  respect  to  the  angle  4>  ,  the  angle  to  the  cell  center 
<$^>  is-  found  by  simply  dividing  the  number  of  cells  in 
ring  i  into  360°.  This  will  provide  an  equal  area  in  each 
cell  of  a  particular  ring. 

pi 

Equation  (4)  can  now  be  used  for  finding  and 

<p  •  > 

Eq  (6)  will  find  ggp  .  The  following  values  for  the 
variables  of  the  ten  cell  model  can  be  found  (Table  XXII). 
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TABLE  XXII 

Ten  Cell  Model  Values 

Ring 

# 

Cell 

Cells  in 
Ring  i 

Cells 

Inside 

Pj/CEP 

<Pi>/CEP 

<4  > 

1 

1 

1 

1 

1 

0.39 

0 

N/A 

2 

2 

4 

S 

1.00 

0.7119 

45° 

2 

3 

4 

S 

1.00 

0.7119 

135° 

2 

4 

4 

5 

1.00 

0.7119 

225° 

2 

5 

4 

5 

1.00 

0.7119 

315° 

3 

6 

5 

10 

00 

1.507 

0° 

3 

7 

5 

10 

00 

1.507 

72® 

3 

8 

5 

10 

00 

1.507 

144® 

3 

9 

5 

10 

00 

1.507 

216° 

3 

10 

5 

10 

00 

1.507 

288° 

It  was  previously  shown  that: 

nt 

PK(X)  -  Jr-  I  PD(I)  . 

nT  i«i 

To  find  PK(X)  ,  PD(I)  must  be  determined.  One  well  accepted 
approach  is  to  define  PD  (I)  as  a  cumulative  log-normal 
function: 
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PD(I) 


I 

/ 

0 


■«*(- 


*nI'-AnI 


T 


.  5  ^  2 


dl 


Here  3  is  the  slope  of  the  intensity  versus  PD  curve 
(plotted  as  a  straight  line  on  log- log  paper)  and  ini  5 
is  the  intensity  at  PD  (I)  *  0.5,  now  called  a  .  Letting 
Z  *  (^~g— ' (therefore  dZ  *  dl)  we  get: 

PD  (I)  -  /  -L-  e'HZl  dZ 

-oo  /2ir 


or  PD  (I)  equals  the  standardized  normal  probability  density 
function  evaluated  from  -«  to  Z  . 

By  defining  sure-kill  (Isk)  and  sure-safe  (Iss) 
probabilities  at  98%  kill  and  2%  kill,  respectively,  the 
area  under  the  normal  curve  is  0.98  and  0.02,  respectively, 
and  the  upper  limits  of  integration  are  just  equal  to  the 
Z  value  in  a  normal  table  corresponding  to  these  areas. 
Therefore: 


*nIsk'a 

g 


+2.054 


and 


fcnl  -a 
ss 


-2.054  . 


Solving  these  equations  simultaneously,  we  get 


a  -  Jta((Isk)(Iss)}  and  B  -  TT2TT5IT)  ' 
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where  I  is  the  sure-safe  intensity  and  IeV.  is  the 


sk 


sure-kill  intensity. 

Now  the  values  of  PK(X)  can  be  computed  as  functions 
of  intensity.  The  required  inputs  are  the  yield  and  CEP  of 
the  attacking  weapons,  the  sure-safe  and  sure-kill  intensities 
(for  the  weapon  effects  of  interest) ,  and  the  distance  of 
the  target  from  the  DGZ.  Given  these  inputs,  the  ten  cell 
model  can  determine  the  probability  of  kill  of  a  weapon 
against  a  target.  In  this  study,  the  weapon  effects  used 
are  neutrons  for  determining  interceptor  PK,  and  overpressure 
and  cratering  for  determining  RV  PK. 
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Appendix  C 


Probability  of  Kill  Routine  for  Overpressure 

To  calculate  the  probability  of  kill  (PK)  of  the  shelters 
due  to  overpressure,  a  procedure  based  on  the  ten  cell  model 
is  used.  The  ten  cell  model  (Appendix  B)  requires  the  sure- 
safe  and  sure-kill  overpressures  of  the  shelter.  The  sure- 
safe  overpressure  is  the  overpressure  at  which  survival  is 
expected  98  percent  of  the  time,  and  the  sure-kill  overpressure 
is  the  overpressure  at  which  destruction  is  expected  98  per¬ 
cent  of  the  time  (Ref  5). 

The  overpressure  caused  by  a  surface  burst  which  creates 
a  crater  is  similar  to  that  of  a  "free-air"  burst  with  one 
and  one-half  times  the  yield  (Ref  21).  The  graph  of  peak 
overpressure  (psi)  versus  distance  from  the  burst  in  feet  or 
meters  for  a  "free-air"  burst  of  a  one  kiloton  weapon  at  sea 
level  is  shown  in  Figure  13  (Ref  3) .  An  equation  approximating 
this  graph  is  as  follows: 

PSI  -  eC . 19 (AnSR) 2  -  1 . 5 (tnSR)  -  .1] 

where 

SR  ■  scaled  distance  from  burst  in  kilometers 

/Atmospheric  Pressure  at  Burst  Altitude  \  ^ 

SR  -  (Actual  Distance)  V  Atmospheric  Pressure  at  Sea  Level  J 

(1.5  Yield) 1/3 

(Ref  25) 
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Table  XXIII  checks  the  difference  between  the  values  read 


from  the  graph  and  the  values  provided  by  the  equation.  The 

equation  provides  an  approximation,  which  on  the  average  is 

/ 

within  S.l  percent  of  the  graph. 


TABLE  XXIII 

Graph  versus  Equation 

Slant  Range 

PSI 

PSI 

(feet) 

(graph) 

(equation) 

Difference 

30 

71,000 

68,151 

2849 

150 

550 

565 

15 

400 

46 

42 

4 

1000 

7 

7 

0 

5000 

.5 

.55 

.05 

The  equation  given  above  provides  the  scaled  overpressure  for 
a  one  kiloton  weapon  at  sea  level.  The  actual  overpressure 
(PSIA)  is  found  using  the  following  equation: 

pct  *  PcT  Atmospheric  Pressure  at  Burst  Altitude\ 
l  Atmospheric  Pressure  at  £ea  Level  j 

Thus,  the  overpressure  of  a  surface  burst  can  be  found  using 
the  following  equation: 
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T  m  [.19(JtaSR)2-1.5(JinSR)-.lII  Atmospheric  Pressure  at  Burst  Altitude 
ir=>iA.  *  e  Atmospheric  Pressure  at  Sea  Level 


where 


SR  =  scaled  distance  from  burst  in  kilometers. 


Appendix  D 


Probability  of  Kill  Routine  for  Neutron  Fluence 

To  calculate  the  probability  of  kill  (PK)  of  the  RVs  due 
to  neutron  fluence,  the  ten  cell  model  is  used  (Appendix  B) . 

The  sure-safe  and  sure-kill  neutron  fluence  levels  of  an  RV 
are  assumed  to  be  1013  and  1017  neutrons  per  square  centimeter 
(N/CM2),  respectively  (Ref  22).  Although  the  actual  fluence 
levels  depend  on  the  RV  design,  these  unclassified  values 
chosen  appear  to  be  reasonable  (Ref  7). 

To  calculate  the  actual  neutron  fluence  at  various  ranges, 
the  number  of  neutrons  produced  per  kiloton  yield  of  the 
terminal  interceptor  must  be  known.  Since  a  thermonuclear 
device  produces  more  neutrons  per  kiloton  than  a  pure  fission 
device,  the  terminal  interceptor  warhead  is  assumed  to  be  a 
50-50  fission-fusion  device.  The  neutrons  per  kiloton  yield 
for  a  50-50  thermonuclear  device  is  approximately  3.16  x  1023 
neutrons  per  kiloton  (Ref  4). 

The  4ttR2  neutron  fluence  versus  mass  integral  (MI) 
where  R  is  the  distance  from  burst  point  to  target  is  shown 
in  Figure  14.  The  mass  integral  (MI)  is  the  amount  of  air 
that  a  neutron  must  traverse  in  traveling  from  the  burst  point 
to  the  target.  An  equation  which  fits  this  graph  is  given 
by: 
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PI*R 


G/CM**2 

Figure  14.  4ttR2  Neutron  Fluence  for  Fission 

and  Thermonuclear  Sources  (Ref  13:55) 
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4irR2  Fluence 


exp  [-6.775  +  .5269x10' 2 (MI)  -  .54364x10' 5 (MI)2 


-  .21468x10' 3 (MI)3/2  -  3.8214 (MI)1/2 

+  10.875(MI)1/3  -  1.3975 (to (MI) )3  (Ref  13:51-55). 

It  is  assumed  that  the  actual  terminal  intercept  will  occur 
at  approximately  20,000  feet  altitude.  Since  the  terminal 
interceptor's  probability  of  kill  is  determined  when  it 
passes  through  the  RV  altitude,  a  homogeneous  atmosphere  with 
a  density  (p)  of  0.65312  kilograms  per  meter  cubed  (KG/M3) 
was  assumed.  This  value  is  used  to  calculate  the  MI  in  a 
homogeneous  atmosphere  as  follows: 

MI  =  pR 

where  R  is  the  distance  from  the  target  to  the  burst  point. 
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Appendix  E 


Analytical  Attack  Probability  of  Kill 
for  Exoatmospheric  Defense  Model 


where 

A  ■  Upper  layer  of  defense  probability  of  kill 

B  ■  Number  of  exoatmospheric  interceptors 

C  ■  Number  of  attacking  RVs 

D  ■  Exoatmospheric  interceptor  probability 
of  kill 

D  -  1  -  exp(-.  694  (WR/CE  ’) 2) 
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where 


A  *  Probability  that  the  MX  survives  the  1st  RV 
targeted  at  it 

B  ■  Probability  that  the  MX  survives  the  2nd  RV 
targeted  at  it 

C  «  Probability  that  the  MX  is  destroyed  by  the  1st 
RV  targeted  at  it 

D  *  Probability  that  the  MX  is  destroyed  by  the  2nd 
RV  targeted  at  it 

E  *  Probability  of  destroying  the  1st  RV  targeted 
at  the  MX 

F  *  Probability  of  the  1st  RV  targeted  at  the  MX 
actually  being  targeted  at  the  MX  (#  of  RVs/# 
of  shelters  or  1,  whichever  is  smaller) 

G  *  Attacking  RV  probability  of  destroying  an  MX 
shelter  (0.5986  for  a  1000  XT  RV  with  a  CEP 
of  0.2  NM) 

H  ■  Probability  of  destroying  the  2nd  RV  targeted 
at  the  MX 

I  ■  Probability  of  the  2nd  RV  targeted  at  the  MX 
actually  being  targeted  at  the  MX  ((#  of  RVs- 
#  of  shelters)/#  of  shelters  or  J0,  whichever 
is  larger) 

J  *  Probability  of  the  terminal  defense  unit  (TDU) 
being  active  for  the  1st  RV  targeted  at  the  MX 

K  ■  Probability  that  the  TDU  has  interceptors  available 
for  the  1st  RV  targeted  at  the  MX 

L  -  Defensive  interceptor  probability  of  kill 

(0.599  for  a  5  KT  interceptor  with  a  CEP  of 
600  feet) 

M  *  Probability  of  the  TDU  being  active  for  the 
2nd  RV  targeted  at  the  MX 

N  »  Probability  that  the  TDU  has  interceptors 

available  for  the  2nd  RV  targeted  at  the  MX 

P  *  Probability  that  the  TDU  does  not  have  inter¬ 
ceptors  available  for  the  2nd  RV  targeted  at 
the  MX 
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*  Probability  of  2  RVs  being  targeted  at  the  TDU 
( (#  of  RVs  -  #  of  shelters)/#  of  shelters  or 
0,  whichever  is  larger) 

*  Probability  of  2  RVs  being  targeted  at  the  MX 
shelter  ((#  of  RVs  -  #  of  shelters)/#  of  shelters 
or  0,  whichever  is  larger) 

■  Probability  of  having  no  opportunity  to  intercept 
the  2nd  RV  targeted  at  the  MX 

■  Probability  that  the  TDU  is  dead  before  the 
1st  RV  targeted  at  the  MX  can  be  intercepted 

*  Probability  that  the  TDU  is  dead  before  the 
2nd  RV  targeted  at  the  MX  can  be  intercepted 

■  Probability  that  the  TDU  is  destoryed  by  both 
RVs  targeted  at  it 

■  Probability  that  the  TDU  survives  the  1st  RV 
targeted  at  it 

«  Probability  that  the  TDU  survives  the  2nd 
RV  targeted  at  it 

■  Probability  that  the  TDU  is  destroyed  by  the 
1st  RV  targeted  at  it 

»  Probability  of  destroying  the  1st  RV  targeted 
at  the  TDU 

*  Probability  of  the  1st  RV  targeted  at  the  TDU 
actually  being  targeted  at  the  TDU  (#  of  RVs/ 

#  of  shelters  or  1,  whichever  is  smaller) 

■  Attacking  RV  probability  of  destroying  a  TDU 
shelter  (0.5852  for  a  1000  KT  RV  with  a  CEP 
of  0.2  NM) 

■  Probability  that  the  TDU  is  destoryed  by  the 
2nd  RV  targeted  at  it 

*  Probability  of  destroying  the  2nd  RV  targeted 
at  the  TDU 

»  Probability  of  the  2nd  RV  targeted  at  the  TDU 
actually  being  targeted  at  the  TDU  ((#  of 
RVs  •  #  of  shelters)/#  of  shelters  or  0,  which¬ 
ever  is  larger) 


DD  ■  Probability  of  the  TDU  being  active  for  the  1st 
RV  targeted  at  the  TDU 

EE  *  Probability  that  the  TDU  has  interceptors 

available  for  the  1st  RV  targeted  at  the  TDU 

FF  *  Probability  of  the  TDU  being  active  for  the 
2nd  RV  targeted  at  the  TDU 

GG  *  Probability  that  the  TDU  has  interceptors 

available  for  the  2nd  RV  targeted  at  the  TDU 

HH  *  Probability  that  the  TDU  does  not  have  inter¬ 
ceptors  available  for  the  2nd  RV  targeted  at 
the  TDU 

II  *  Probability  of  having  no  opportunity  to  inter¬ 
cept  the  2nd  RV  targeted  at  the  TDU 

JJ  *  Probability  that  the  TDU  is  destoryed  by  the 
1st  RV  targeted  at  the  TDU 
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Appendix  G 

Computer  Model  Listings 
and 

Q-GERT  Networks 
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41 


Computer  Model  Listin 


Q-GERT  Listing  (Layered  Defense  Model) 

GEN*SH£RET  vLAYERl *11*6*1 98 1*5*0*0*20. *600*S*0. >4* 

S0U*1»0*1*A»M* 

VAS*  1 » 1 »IM» 1*2* CO* 10»4*C0» 12*  4C0UNTER » SHEL-2 » SHELTERS  * 

ACT *1*1* NO *1* 1 /GENERATE*  (9>A1 .LT .16*  GENERATE  RV  ARIV  * 

PAR  *  1  *  » 083  *  0 .  *  *  . 033  *  94  RVS  ARRIVE  NORMALLY .  MEAN-5  3EC*SD«2  SEC 
ACT*1*3*C0*0«  *2/RVARIV* <9)A1 .LE.16*  RVS  ARRIVE  TO  SYSTEM  * 

REG* 3* 1 » 1 »F*  DETERMINE  IF  RVS  ARE  TWICE  *  OF  SHELTERS 
ACT*3*41 *C0* 1 » » 3/TWOPER* ( 9> A4.EQ.8#  RVS  TWICE  *  OF  SHELTERS  * 
ACT * 3* 30* CO* 0*  *4/LTTW0PER*  <9>A4.GT»8*  RVS  LESS  THAN  TWICE  « 

REG»41 *  1  *  1  »P*  TWICE  *  OF  RVS  ASSIGNED  TO  SHELTERS 
ACT *41*42*C0»0«  *5/MPS* (8) .8333*  ASSIGN  RVS  TO  MPS  * 

ACT*41*43*C0*0.*5/DU*<8>. 08335*  ASSIGN  RVS  TO  DU  * 

ACT *41 *44*C0*0»  *5/MX* (8>  *08335*  ASSIGN  RVS  TO  MX  * 

QUE * 42/MPSQUE * 0 * 0 * D * F *41*  MPS  QUEUE  # 

QUE»43/DUQUE*0*0*D*F *41*  DU  QUEUE 
QUE»44/MXQUE*0*0*D»F»41*  MX  QUEUE 
VAS*42*3*C0* 1*  ASSIGN  MPS  ID  * 

VA3 *43*3* CO *2*  ASSIGN  DU  ID  * 

VAS  *  44  *  3  *  CO  *  3*  ASSIGN  MX  ID  # 

ACT*42* 19*C0* 10 .  *7/MPSDELAY *  12*  SEND  MPS  RVS  TO  MAIN  RV  QUEUE(RV-4 

ACT *43* 19* CO *10.  * 8/DUDELAY»2*  SEND  DU  RVS  TO  MAIN  RV  QUEUE 

ACT *44* 19* CO* 10* » 9/MXDELAY»2*  SEND  MX  RVS  TO  MAIN  RV  QUEUE 

QUE »30*0**D*F«  LESS  THAN  TWICE  THE  *  OF  RVS  QUEUE 

ACT  * 30 *31* CO *1*  22/CAPCK  *  24*  CHECK  MPS/DU/MX  QUE  STATUS 

REG»31 *  1  *  1 *P*  LESS  THAN  TWICE  THE  ♦  OF  RVS  ASSIGNED  TO  SHELTERS 

ACT *31 *32*C0*0.  *6/MPS* (8) *8333*  ASSIGN  RVS  TO  MPS  • 

ACT»31 »33*C0*0.  * A/DU* (8) *08335*  ASSIGN  RVS  TO  DU  « 

ACT *31*34*C0*0»  *6/MX*  <8) *08335*  ASSIGN  RVS  TO  MX  * 

QUE  *  32/MPSQUE  *0*0»D*F*31«  MPS  QUEUE  « 

QUE  *  33/DUQUE *0»0*D*F*31«  DU  QUEUE 

QUE*34/MXQUE»0»0*D*F *31*  MX  QUEUE 

ACT *32* 19* CO* 10*  * 10/MPSDELAY *14*  SEND  MPS  RVS  TO  MAIN  RV  QUEUE 
ACT *33* 19* CO *10*  *  1 1/DUDELAY* 2*  SEND  DU  RVS  TO  MAIN  RV  QUEUE 

ACT*34* 19*C0* 10«  * 12/MXDELAY  »2*  SEND  MX  RVS  TO  MAIN  RV  QUEUE 

VAS*32*3*C0» 1 *2*US*2«  ASSIGN  MPS  ID  • 

VAS  *  33  *  3  *  CO  *  2  *  2  »  US  *  2*  ASSIGN  DU  ID  * 

VAS  *  34  *  3  *  CO  »  3  *  2  *  US  *  2*  ASSIGN  MX  ID  * 

QUE* 19/RVQUE*0*  *D*F«  MAIN  RV  QUEUE 
VAS»19*4»C0» 1*  ASSIGN  ALIVE  (1>  CODE 
ACT  * 19*20*C0*O>  * 13/RVDELAY*  RV  DELAY 
REG*20*1*1*F«  DETERMINE  RVS  TO  SEND  TO  PROBE 

ACT*20*21 *C0*0»  *14/T0PR0BEr (9)A1 <LE»5«  SEND  RVS  TO  PROBE  QUE  * 
ACT *20*51* CO *»5* 15/T0L0ADS* (9)A1 .GT.5*  EXCESS  RVS  FOR  STAT  COL  * 
QUE*51/RVACTIVE*0*  *D»F* 

QUE*21/PR0BEQUE*0» »D*F*  *  10. * *5«  PROBE  QUEUE 

ACT*21 *23*US*3» 17/INTERCEPT0R»5*  INTERCEPTOR/RV  INTERACTION  * 
PAR*2*0«  *  *  *1*6988*9*  XIY  DISTRIBUTION 
STA*23/P1STATS*1*1*D*I*10.*.5«  PROBE  1  STATS 
ACT  *  23  *  24  *  CO  *  0  *  18/DELAY* 
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QUE  *  24/CKSTATUS » 0  * » D  »  F* 

ACT *24*50*US* 1  * 19/CKALIVE*5#  SEND  ACTIVE  RVS  TO  NODE  51  * 

STA*50/DEADSTAT*1*1»D*I»10. * .5*  COLLECTS  STATS  ON  KILLED  RVS 
ACT *  51  *  53  *  CO  *  3 . * 20/LDDELAY >24*  DELAY  UNTIL  LOADS  INTERACTION 
QUE  » 53/LDSQUE *  0  *  *  D » F*  LOADS  QUEUE 
ACT*53»54»US»4»21/LDSSERV*24*  LOADS  INTERACTION 
REQ  *  54  *  1 » 1 »  D* 

ACT»54»55»US»5*23/MXSTATUS*24*  CK  STATUS  OF  RVS  &  HX 
STA»55/LDSTATS» 1* 1*D*I*13. » .5*  LOADS  STATISTICS 
FIN* 


Q-GHVT  Listing  C&idoatmospheric  Defense  Model  -  One  TDU) 


GEN 

SOU 

VAS 

ACT 

PAR 

ACT 

REG 

ACT 

ACT 

REG 

ACT 

ACT 

ACT 

QUE 

QUE 

QUE 

VAS 

VAS 

VAS 

ACT 

ACT 

ACT 

QUE 

ACT 

REG 

ACT 

ACT 

ACT 

QUE 

QUE 

QUE 

ACT 

ACT 

ACT 

VAS 

VAS 


SHERET»LOADS*ll» 16* 1981 » 1 » 0*0*20 . *600»S»0.»4* 

1 » 0  * 1 * A  * M* 

1*1*IN*1 *2tCQt 10* 4* CO* 12*  *COUNTER * SHEL-2 * SHELTERS 

1  *  1 *NG* 1 *  1 /GENERATE* <9>A1.LT.16*  GENERATE  RV  ARIV 
it .083*0. * * .033*6*  RVS  ARRIVE  NORMALLY »  MEAN=5  SEC»SD=2  SEC 
l*3»C0f0«  »2/RVARIVr (?>A1«LE»16*  RVS  ARRIVE  TO  SYSTEM 
3* 1 r 1 »F*  DETERMINE  IF  RVS  ARE  TWICE  *  OF  SHELERS 
3*41 rCOr 1 »  » 3/TUOPER» (9) A4« EQ»8«  RVS  TWICE  *  OF  SHELTERS 
3»30*C0»0»* 4/LTTWOPER » ( 9 ) A4 . GT . 8*  RVS  LESS  THAN  TWICE 
41'1'ltP*  TWICE  *  OF  RVS  ASSIGNED  TO  SHELERS 
41* 42»C0* O*  * 5/MPS* <8>  *8333*  ASSIGN  RVS  TO  MPS 
41 *43»C0*0.*5/DU»<8). 08335*  ASSIGN  RVS  TO  DU 
41*44*C0f0«» 5/MX » < 8 ) ♦ 08335*  ASSIGN  RVS  TO  MX 
42/MPSQUE »0*0fDfF»41*  MPS  QUEUE 
43/DUQUE»0* 0»D»F»41«  DU  QUEUE 
44/MXQUE »0»0fD»F»41*  MX  QUEUE 
42»3'C0*1*  ASSIGN  MPS  ID  * 

43*3* CO* 2*  ASSIGN  DU  ID  * 

44 *3* CO* 3*  ASSIGN  MX  ID  * 

42* 1?»C0» 10. »7/MPSDELAY *12*  SEND  MPS  RVS  TO  MAIN  RV  QUEUE(RV-4) 
43* 19* CO *10. *8/DUDELAY *2*  SEND  DU  RVS  TO  MAIN  RV  QUEUE 
44* 19* CO* 10. »9/MXDELAY»2*  SEND  MX  RVS  TO  MAIN  RV  QUEUE 
30*0* »D»F*  LESS  THAN  TWICE  THE  *  OF  RVS  QUEUE 
30 *31* CO *1.* 22/C APCK * 24*  CHECK  MPS/DU/MX  QUE  STATUS 
31  *  1  *  1 *P«  LESS  THAN  TWICE  THE  *  OF  RVS  ASSIGNED  TO  SHELTERS 
31 *32*C0*0. »6/MPS* < 8) .8333*  ASSIGN  RVS  TO  MPS  * 

31* 33* CO *0.*  6/DU  *  <  S  > • 08335*  ASSIGN  RVS  TO  DU  * 

31* 34 * CO *0.*  4/MX  * ( 8  > • 08335*  ASSIGN  RVS  TO  MX  * 

32/MPSQUE *0*0*D»F*31*  MPS  QUEUE  * 

33/DUQUE*0*0*D*F*31»  DU  QUEUE 
34/MXQUE*0*0*D*F*31«  MX  QUEUE 

32* 19*C0* 10. * 10/MPSDELAY *  14*  SEND  MPS  RVS  TO  MAIN  RV  QUEUE 
33* 19*C0* 10. * 11/DUDELAY * 2*  SEND  DU  RVS  TO  MAIN  RV  QUEUE 
34*  19»C0*  10.  *  12/MXDELAY * 2*  SEND  MX  RVS  TO  AN  RV  QUEUE 
32 *3* CO *1*2*  US  *  2*  ASSIGN  MPS  ID  * 

33  *  3  *  CO  *  2  *  2  *  US  *  2*  ASSIGN  DU  ID  • 
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VAS 1 34 » 3 » CO » 3 » 2 1 US » 2*  ASSIGN  MX  ID  * 

QUE  f 19/RVQUEf  Of f  Df  F*  MAIN  RV  QUEUE 

VASf 19f4fC0fl*  ASSIGN  ALIVE  <1)  CODE 

ACTfl9f53fC0f3.f20/T0L0ADSf24*  SEND  RVS  TO  LOADS 

QUE»53/LDSQUErO»  >D»F*  LOADS  QUEUE 

ACT f53f 54fUSf 4f 21/LDSSERVf 24*  LOADS  INTERACTION 

REGf 54f 1 f 1 f  D* 

ACT»54»55»US»5»23/MXSTATUS»24*  CK  STATUS  OF  RVS  &  MX 
STAfS5/LDSTATS»lf l»Df Ifl3. » .5*  LOADS  STATISTICS 
FIN* 


j-GBIT  Listing  (Endoatmospheric  Defense  Model  -  Two  TDUs) 

3EN»SHERET»LDS2DU»11»20» 1981 f2f Of Of  35. f 1200fSf0. »5* 

»0U»1»0»1 »A»M* 

/ASflf It IN» l»2»C0f9»4»C0f 12*  #COUNTER * SHEL-3 1 SHELTERS  * 

*CT flflfNOflf 1/GENERATE » (9)A1.LT*16#  GENERATE  RV  ARIV  * 

^AR » 1 r » 083 f  0 » f  f  « 033 » 3*  RVS  ARRIVE  NORMALLY ♦  MEAN=5  SEC»SD=2  SEC 
ICT r 1 »3 rCO »0 . f 2/RVARIVf ( 9) A1 . LE » 16*  RVS  ARRIVE  TO  SYSTEM  * 

tEGr3f 1 t 1 f F*  DETERMINE  IF  RVS  ARE  TWICE  *  OF  SHELERS 
*CT r3»41 »C0» 1 . f 3/TWOPER? (9) A4.EQ.8*  RVS  TWICE  *  OF  SHELTERS  * 

^CT » 3 1 30 » CO » 0 . » 4/LTTWOPER f ( 9 ) A4 ♦ GT  *  8*  RVS  LESS  THAN  TWICE  * 

MEG  f  41 » 1 » 1 »P*  TWICE  *  OF  RVS  ASSIGNED  TO  SHELERS 
)CTf 41 f 42f COf 0. f 5/MPSf  <8) .75*  ASSIGN  RVS  TO  MPS  (12  SHEL) 

ICT r41»43»C0»0.» 5/DU » ( 8 >  « 0833333334*  ASSIGN  RVS  TO  DU  1  (12  SHEL) 

\CT » 41 » 44 fC0f0.f5/MXf (8). 0833333333*  ASSIGN  RVS  TO  MX  (12  SHEL) 

ICT * 41 ' 45 f CO r0.f5/DU2' (8) *0833333333*  ASSIGN  RVS  TO  DU2  (12  SHEL) 

1UE f 42/MPSQUE »0»0»D»F»41*  MPS  QUEUE  * 

)UEf 43/DUQUE»0f Of DoF  f  41*  DU  1  QUEUE 
lUEt 44/MXQUEf Of Of D» Ff 41*  MX  QUEUE 
)UEf45/DU2QUEf0f0fDfFf41*  DU2  QUEUE 
)ASf 42f3f COf 1*  ASSIGN  MPS  ID  * 

FASf 43f 3f COf 2*  ASSIGN  DU  1  ID  * 

^ASf44f3fC0f3*  ASSIGN  MX  ID  * 

)ASf 45f 3f COf 4*  ASSIGN  DU2  ID  * 

kCT s*42f  19f  COf  10*  f  7/MPSDELAY f  10*  SEND  MPS  RVS  TO  MAIN  RV  QUEUE (RV-6) 
kCT f43fl9fC0fl0»  f 8/DUDELAY 1 2*  SEND  DU  RVS  TO  MAIN  RV  QUEUE 
)CT f44fl9fC0fl0. f 9/MXDELAYf 2*  SEND  MX  RVS  TO  MAIN  RV  QUEUE 
*CTf45f 19fC0fl0.  99/DUDELAY 1 2* 

)UEf 30t Of  f Df F*  LESS  THAN  TWICE  THE  *  OF  RVS  QUEUE 
kCTf 30f 31 f COf 1  *  f 22/CAPCKf 24*  CHECK  MPS/DU/MX  QUE  STATUS 
(EGf 31 f 1 f 1 f P«  LESS  THAN  TWICE  THE  *  OF  RVS  ASSIGNED  TO  SHELTERS 
ICT f31f32f COf 0*  f 6/MPSf (8) >75*  ASSIGN  RVS  TO  MPS  (12  SHEL) 

ICTf 31 f 33f COf 0*  f 6/DUf (8) *0833333334*  ASSIGN  RVS  TO  DU  1  (12  SHE 

kCT f 31 f 34f COf 0*  f 6/MXf (8) *0833333333*  ASSIGN  RVS  TO  MX  (12  SHEL) 

*CTf 31 f 35 f COf O.f A/DU f (8). 0833333333*  ASSIGN  RVS  TO  DU2  (12  SHEL) 
)UEf32/MPSQUEf0f0fDfFf31*  MPSQUE  * 

)UEf33/DUQUEf0f0fDfFf31*  DU  QUEUE 
lUEf 34/MXQUEf Of Of Df F f 31*  MX  QUEUE 


(12  SHEL) 
(12  SHEL) 
(12  SHEL) 

* 
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QUE * 3S/DU2QUE ,0»0*D,F,31*  0U2  QUEUE 

ACT *32*19*C0*10«* 10/MPSDELAY ,13*  SEND  MPS  RVS  TO  HAIN  RV  QUEUE(RV-3> 
ACT *33 *19* CO ,10*, 11/DUDELAY * 2*  SEND  DU  RVS  TO  MAIN  RV  QUEUE 
ACT*34, 19, CO, 10.  * 12/MXDELAY,2*  SEND  MX  RVS  TO  MAIN  RV  QUEUE 
ACT*35* 19*C0* 10*  ,98/DUDELAY* 2*  SEND  DU2  RVS  TO  MAIN  RV  QUEUE 
VAS  *  32  »  3  »  CO  *  1  *  2 » US  *  24c  ASSIGN  MPS  ID  * 

VAS  *33*3*  CO  *  2  *  2  *  US  *  2*  ASSIGN  DU  1  ID  * 

VAS  *  34  *  3  *  CO  *3*2* US *2*  ASSIGN  MX  ID  * 

VAS  *  35  *  3  *  CO  *  4  *  2  *  US  *  2*  ASSIGN  DU2  ID  # 

QUE  *  1 9/RVQUE *0  **D*F»  MAIN  RV  QUEUE 

VAS  *  1 9  *  4  *  CO  *  1 4  ASSIGN  ALIVE  (1)  CODE 

ACT  * 19*53*C0* 1  *  *20/T0L0ADS  *  24#  SEND  RVS  TO  LOADS 

REG  *  53  *  1  *  1  *  D* 

ACT*53»54»C0*0.,31/LDSSERV,24*  LOADS  INTERACTION 
REB*54*1*1*D« 

VAS  *  54  *  5  *  US  *  44 

ACT *54»57*US*5*32/MXSTATUS * 24#  CK  STATUS  OF  RVS  &  MX 
STA,57/LDSTATS»1*1*D*I*13. *.5*  LOADS  STATISTICS 
REG *55* 1  *  1  *D* 

ACT  *  55  *  56  *  CO  *  0  *  *  33/DU2SERV  *  244  RV/DU2  INTERACTION 
REG *56* 1  *  1  ,D* 

VAS, 56, 5, US, 6* 

ACT »56»58»US*5»34/MXSTATUS*24#  DETERMINE  MX  STATUS  AFTER  ATTACK 

STA*58/DU2STAT*1*1*D*I*13. * .5*  DU2  STATS 

FIN* 


SUBROUTINE  US<ISN»DTIM) 

COMMON/USER/MISS * X  »  Y  *  ZZ * DU » DUPK  , MXPK »  J I »  L  *  PK 
COMMON/QVAR/NDE * NFTBU < 1 00 ) » NREL ( 1 00  > » NRELP < 1 00  >  r NREL2 < 1 00 ) » 
INRUN  * NRUNS  *  NTC  < 1 00 ) » P ARAM  < 100*4) *TBEG»TNOW 
C**DECLARE* VARIABLES  ** 

INTEGER  DU»JI»L*I*JK 

REAL  MI SS  »  ZZ  »  X » Y » NO ,  DUPK  * HXPK  » PK  » RA » RE » RCPK # NALPHA » NBETA 
REAL  PIR2F»HI»PI»PERKT,SR#ATT(4>»DENSfF»C(?>»Z»CPK 
REAL  SOP t OP t CEP f YIELD* OPALPHA * OPBETA * T » YLD t TR » MXPS 
DATA  JK*CPK/Of 0./ 

Z-0. 

GO  TO  (1*2*3*4*5)» ISN 

1  DTIH-1. 

ATT  < 1 ) “GATRB ( 1 > 

ATT ( 2 ) “GATRB ( 2 ) 

ATT  <  3 ) “GATRB  <  3 ) 

ATT  <  4 ) *G ATRB <  4 ) 

IF(GATRB<4) . GT. .5)  THEN 
CALL  STAGO( 19»51»0.fl »ATT) 

ENDIF 

RETURN 

C**  ENSURES  CORRECT  RU  TARGETING »  *  SHELTERS  >  1/2  *  RVS. 

2  DTIH-O. 

TRBREHST < 22 ) 

IF< <NTC<32) . GE ♦ 1 3 ) .AND. (NTC(34> .EQ.2) )  THEN 
CALL  ST  AGO  <22»33»TRfO»ATT) 

ELSE IF ( <NTC(32) .GE. 13) . AND. <NTC<33) .EQ.2) )THEN 
CALL  ST AGO< 22*34 *TR»0» ATT) 

ELSEIF ( <NTC(32) .GE. 14) . AND. <NTC<34> .EQ. 1 ) )  THEN 
CALL  STA5O  <  22  *  33  »  TR  *  0  *  ATT  > 

ELSEIF ( (NTC(32) .GE.14) .AND. <NTC<33) .EQ.l) )  THEN 
CALL  ST  AGO  <  22 » 34  *  TR  »  0  »  ATT ) 

ELSEIF <NTC< 32 ) .EQ. 14)  THEN 
RN-DRAND(9) 

IF<RN.LE. .5)  THEN 

CALL  ST AGO ( 22  *  33  *  TR  »  0  *  ATT ) 

ELSE 

CALL  ST AGO ( 22  *  34  *  TR  *  0  *  ATT ) 

ENDIF 

ENDIF 

RETURN 

C**  SEND  RU  TO  INTERCEPTORS  ** 

C**  FIRE  INTERCEPTOR  AT  RV  I  DETERMINE  HIT/HISS  ** 

3  X-N0<2) 

Y-N0<2) 

ZZ«<X*X+Y»Y)**.5 

HISS-ZZ-1.5 

HISS-ABS(MISS) 
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C**  IS  RV  WITHIN  INTERCEPTOR  WEAPON  RADIUS?  (KILL)  ** 
IF(MISS .LE .1.5)  THEN 
CALL  PATRB(0»4) 

ENDIF 

C**  30  SEC  FROM  DETECTION  TILL  INTECEPTOR/RU  ENCOUNTER  ** 
DTIM-.5 
RETURN 

C**  SEND  RUS  NOT  KILLED  TO  LOWER  LAYER  (LOADS)  ** 

C**  IS  RV  TARGETED  AT  MPS?  ** 

4  IF  (GATRB(3) . EQ. 1 .  )  THEN 
DTIM-O. 

RETURN 

ENDIF 

C**  CHECK  STATUS  OF  DU  ** 

IF(DU.GT..5)  THEN 
DTIM-. 08333 
GO  TO  100 
ENDIF 

IF  (JI.GE.2)  THEN 

C**  IS  THIRD  RV  TARGETED  AT  DU?  ** 

IF  ( GATRB ( 3 ) ♦ EQ . 2 . )  THEN 
DTIM-. 08333 
GO  TO  100 
ENDIF 
ENDIF 

C**  ALL  INTERCEPTORS  LAUNCHED  ** 

DTIM-. 08333 
IF ( JI .GE.3)  THEN 
GO  TO  100 
ENDIF 
JI*JI+1 

C**  INTERCEPTOR  NEUTRON  PK  ** 

CttGLOSSARY** 

C*  MI-MASS  INTEGRAL  (G/CM2) 

C*  PIR2F»4*PI*R*  2  FLUENCE 

C*  PERKT-NEUTRONS  PER  KILOTON 

C*  INTALT-ALTITUDE  OF  INTERCEPT  (KM) 

C*  SR-DI STANCE  BETWEN  BURST  AND  TARGET  (KM) 

C*  DENS-AIR  DENSITY  AT  INTERCEPT  ALTITUDE  (KG/M3) 

C*  YIELD-THERMONUCLEAR  YIELD  (KT) 

C*  NALPHA-NEUTRON  PK  FUNCTION  ALPHA 

C*  NBET A-NEUTRON  PK  FUNCTION  BETA 

C*  CEP-INTERCEPTOR  CIRCULAR  ERROR  PROBABLE 
C*  F-NEUTRON  FLUENCE  ( NEUTR0NS/CM2 ) 

C* 

C**  INPUT  DATA  t* 

CEP-600 
PI-3.1415927 
PERKT-3 . 16E23 
DENS-. 653 12 


"WPl^iiUPipp 


YIELD-5. 

DATA  C/- . 6775E1 > . 5269E-2  > - . 54364E-5  r - . 21 468E-3 *  - . 3821 4E1 1 
i. 10875E2»- . 13975E1/ 

NALPHA-34 • 539 
NBETA-2.242 

C**  CALCULATE  SLANT  RANGE  USING  10  CELL  MODEL  ** 

DO  200  I-l t 10 

CALL  TENCELL ( SR  r  CEP  r I > 

C**  CALCULATE  MASS  INTEGAL  ** 

IF(SR.LT . . 1 )  THEN 
SR*. 001 
ENDIF 

C**  CHANGE  SR  TO  KILOMETERS 
SR-SR*. 0003048 
MI=(DENS)*SR 
MI— MlfclOO. 

C**  COMPUTE  HOMOGENEOUS  AIR  4PIR2  FLUENCE  ** 

PIR2F-EXP<C<1)+C<2)*MI+C(3)*MI**2.+C<4>*MI**1.5+C<5>*MI**.5+ 
*C(A)*MI**(l./3.)+CC7)*<AL0G<MI) )) 

C**  COMPUTE  FLUENCE  ** 

F»<  <PIR2F)/<4.*PI*<SR**2. >  >*(PERKT>*< YIELD) *( IE-10) ) 

C**  Z-NORMAL  TABLE  ENTRY  VALUE  FOR  NEUTRON  FLUENCE  ** 

Z* ( ALOG  <  F ) -NALPHA ) /NBETA 
CALL  PROB(Z.CPK) 

200  CONTINUE 

RN-DRAND (9) 

IF(RN.LE.CPK)  THEN 
CALL  PATRB<0»4) 

CPK-O.  . 

RETURN 

ENDIF 

CPK-O. 

C 

C**  RV  KILL  MECHANISM  ** 

100  DTIM-DTIM+. 08333 

IF ( GATRB (31.EQ.3.)  THEN 
L-L+l 
ENDIF 
C 

C**  RV  OVERPRESSURE  PK  ** 

C**GLOSSARY«* 

C*  OPALPHA-OVERPRESSURE  PK  FUNCTION  ALPHA 

C*  OPBETA-OVERPRESSURE  PK  FUNCTION  BETA 

C*  SOP-SCALED  OVERPRESSURE 

C*  OP-OVERPRESSURE 

C 

C**  INPUT  DATA  ** 

CEP-1215.2 

YLD-1000. 

OPALPHA-6 . 8755 
OPBETA-. 12435 
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C**  CALCULATE  SLANT  RANGE  USING  TEN  CELL  MODEL  ** 

DO  300  I-l»10 

CALL  TENCELL ( SR » CEP  >  I ) 

C**  SCALE  SR  FOR  PRESSURE  AND  YIELD  ** 

SR»SR*< (12. 69/14. 7>**<l,/3.) >/< < 1 . 5*YLD>**< 1 ./3. > > 
C**  CHECK  FOR  CELL  1  ** 

IF  (SR.LT. .1 )  THEN 

S0P-1E5 

ELSE 

C**  CHANGE  SR  TO  KILOMETERS  ** 

SR-SR*. 0003048 

SOP-EXP  < » 19*  < ALOG  <  SR ) ) **2-1 . 5*AL0G  <  SR ) - ♦ 1 ) 

ENDIF 

C**  SCALE  OVERPRESSURE  FOR  ALTITUDE  ** 

OP-SOP* < 12 ♦ 69/14. 7) 

C**  Z-NORMAL  TABLE  ENTRY  VALUE  FOR  RV  OVERPRESSURE  ** 

Z»  <  ALOG ( OP ) -OPALPHA ) /OPBETA 
CALL  PROB<Z»CPK) 

300  CONTINUE 
ROPK-CPK 
CPK-0  » 

C 

C**  RV  CRATERING  PK  ** 

C»*GLOSSARYS* 

C*  RA-APPARENT  CRATER  RADIUS 
C*  RCPK-RV  CRATER  PROB  OF  KILL 
RA— 484.54 
RE-1041.76 
RE-RE* < 2. /3.) 

IF  <  GATRB ( 3 ) . EG « 2 )  THEN 
RCPK-.1704 

ELSEIF  < GATRB ( 3 ) . EQ . 3 )  THEN 

RCPK-. 1971 

ENDIF 

PK-ROPK+ < 1 -ROPK ) *  <  RCPK > 

C**  COMPUTE  CORRECT  PK  FOR  DU  SHELTER  ** 

I F ( GATRB <3).EQ.2)  THEN 

DUPK-PK 

RN-DRAND  <  9 ) 

C**  IS  DU  DESTROYED?  ** 

IF<RN.LE.DUPK>  THEN 
DU-'l 

- - RETURN 

ENDIF 

RETURN 

ENDIF 

C**  COMPUTE  CORRECT  PK  FOR  MX  SHELTER  ** 
IF(GATRB<3>.EQ.3>  THEN 
IF(L.EQ.2)  THEN 
NXPK-l-( 1-MXPK>**2 
ELSEIF(L.EQ.l)  THEN 
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MXPK-PK 

ENDIF 

RETURN 

ENDIF 

C**  IS  HX  DESTROYED?  ** 

5  IF  <NTC<54) »EQ  »NTC<51 > )  THEN 

JI-0 
L-0 

RN-DRAND<9> 

IF<RN.LE.MXPK>  THEN 

PRINT** '  MX  DESTROYED' 

JK-JK+1 

IF  < NRUN » 6E . NRUNS )  THEN 
MXPS-1 ♦ - C  FLOAT  < JK ) /FLOAT  <  NRUNS ) > 

PRINT**  '»>  MX  SURVIVABILITY*  '*MXPS 
ENDIF 
ELSE 

IF < NRUN. GE. NRUNS)  THEN 

MXPS-1 . - ( FLOAT  <  JK ) /FLOAT  <  NRUNS ) ) 

PRINT** '»>  MX  SURVIVABILITY"  '»MXPS 
ENDIF 
ENDIF 
ENDIF 
RETURN 
END 

SUBROUTINE  UI 

COMMON/USER/MISS  *  X  *  Y  *  ZZ  *  DU  *  DUPK  *  MXPK  *  JI *  L  *  PK 
INTEGER  DU*JX*L 

REAL  MISS*X*Y*ZZ*  DUPK  *  MXPK  *  PK  *  CPK  *  Z 
MISS-O. 

X-Y-ZZ-O. 

DU-0 

DUPK-O. 

MXPK-O. 

JI-L-0 

PK-O. 

CPK-0  * 

Z-O. 

RETURN 

END 

SUBROUTINE  TENCELL  <SR*CEP*I> 

C 

COMMON/USER/MISS  *  X  *  Y  *  ZZ  *  DU  *  DUPK  *  MXPK  *  JI *  L  *  PK 
C 

C**GLOSSARY«* 

C*  THETA < I) -REFERENCE  ANGLE  TO  RADII  TO  CELL  CENTERS  MEASURED 

C*  FROM  MAJOR  AXIS  COUNTER  CLOCKWISE 
C*  RHO< D-CELL  RADII  (RHO/CEP) 

C 

C**DECLARE  VARIABLES 
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INTEGER  I 

REAL  THETA < 10) fRH0( 10) fCEPfSR 
C**  DATA  INPUT  *# 

DATA  THETA/0 , 45 f 1 35 f 22S  f315f0f72f144f21<Sf  288/ 

DATA  RH0/0f .7109f  .7109f  .7109 f . 7109 f 1 .509f 1 . 509f 1 . 509f 1 .509f 
*1.509/ 

C**  CALCULATE  INTERCEPTOR  SLANT  RANGE  ** 

SR*RHO  < I ) #CEP 
RETURN 
END 
C 

SUBROUTINE  PROB<ZfCPK> 

COMMON/USER/MISS  f  X  f  Y f ZZ  f  DU  f DUPK  t MXPK  f  JI » L f PK 
C 

CtSGLOSSARYt* 

C*  PK “PROBABILITY  OF  KILL 

C*  CPK-CUM  PROBABILITY  OF  KILL 

C**  DECLARE  VARIABLES  ** 

REAL  PKfCPKfZ 
IF(Z.LT .0. )  THEN 
Z«ABS<Z> 

CM*  CURVE  FIT  FOR  STANDARD  NORMAL  CURVE  ** 

PK* . 5*< 1 . + . 19485*Z+ . 1 15194*Z*Z+ . 000344*Z*Z*Z+ . 019527*Z**4 . ) *#  < -4 . > 
ELSE 

PK-1- . 5*  < 1 .  +  ♦ 1?6854*Z+ . 1 151 94*Z*Z+ . 000344*Z**3 .  + . 01 9527*Z**4 . > 
***(-4.) 

ENDIF 

PK*PK/10. 

CPK*CPK+PK 

RETURN 

END 


USER  Subroutines  (Endoatmospheric  Defense  Model  -  One  TDU) 


SUBROUTINE  U3<ISNfDTIM) 

COMMON/USER/M I SS  f  X  t Y  »  ZZ  f  DU  t DUPK  f  MXPK  »  JI f  L  r  PK 
COMMON/QVAR/NDE  f  NFTBU  < 100  >  f  NREL  < 100 ) f  NRELP  < 100 ) f  NREL2 ( 100 ) f 
INRUN f  NRUNS f  NTC  U  00 ) f  PARAM  < 1 00  f  4 ) f  TBEG  f  TNOW 
INTEGER  DUfJIfLfIfKK 

REAL  MISS f  ZZ  f  X  f  Y  f  NO f  DUPK  f MXPK  f PK  f RA  f RE  f RCPK  f  NALPHA  f  NBET A 
REAL  PIR2FfMX  fPI fPERKT fSRf ATT<4 >  fDENSfF fC<7) fZfCPK 
REAL  SOP f OP f CEP fYIELDf OPALPHA f OPBETA f YLD f TR f MXPS 
DATA  CPKfKK/O.fO/ 

GO  TO  <1f2f3f4f5>fISN 
1  DTZM-1. 

RETURN 

c**  ENSURES  CORRECT  RV  TARGETING.  *  SHELTERS  >  1/2  *  RVS. 
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DTIM-.001 
TR-REMST (22) 

IF < (NTC(32) .GE.13) .AND. <NTC<34> ,EQ.2>  >  THEN 
CALL  STAG0<22t 33»TR»0» ATT ) 

ELSEIF  < <NTC(32) .GE. 13) .AND. <NTC<33> .EQ.2) )THEN 
CALL  ST  AGO  <  22 » 34  »  TR » 0  ?  ATT ) 

ELSE IF ( (NTC<32>  .GE.14)  .AND.  <NTC<34) »EQ»1) >  THEN 
CALL  STAG0<22>33»TR>0»ATT) 

ELSEIF ( ( NTC ( 32 ) .GE.14) .AND. (NTC(33> *EQ*1>>  THEN 
CALL  ST  AGO  <  22  r  34  r  TR » 0 »  ATT  > 

ELSEIF < NTC (32). EG .14)  THEN 
RN-DRAND (6) 

IF(RN*LE. *5)  THEN 

CALL  ST  AGO  <  22  *  33  r  TR  >  0 1  ATT ) 

ELSE 

CALL  ST AGO <  22  f  34  »  TR » 0  »  ATT ) 

ENDIF 
ENDIF 
RETURN 

3  DTIM-1. 

RETURN 

C**  SEND  RV  TO  LOADS  INTERCEPTORS  ** 

C**  IS  RV  TARGETED  AT  MPS  ?  ** 

4  IF  ( GATRB <  3 ) . EQ • 1 . >  THEN 

DTIM-O. 

RETURN 

ENDIF 

C*«  CHECK  STATUS  OF  DU  ** 

IF (DU.GT . .5)  THEN 
DTIM-. 08333 
GO  TO  100 
ENDIF 

IF  <JI.GE.2)  THEN 

C**  IS  THIRD  RV  TARGETED  AT  DUT  ** 

IF  (GATRB<3) .EQ.2. )  THEN 
DTIM-. 08333 
GO  TO  100 
ENDIF 
ENDIF 

C**  ALL  INTERCEPTORS  LAUNCHED  ** 

DTIM-. 08333 
IF ( JI .GE .3)  THEN 
GO  TO  100 
ENDIF 
JI-JI+1 

C**  INTERCEPTOR  NEUTRON  PK  t* 

Ct*GL0S8ARY*t 

C*  HI-MASS  INTEGRAL  (G/CH2) 

C*  PIR2F-4*PI*R*  2  FLUENCE 
C*  PERKT-NEUTRONS  PER  KILOTON 
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C*  INTALT-ALTITUDE  OF  INTERCEPT  (KM) 

C*  SR-DI STANCE  BETUEN  BURST  AND  TARGET  (KM) 

C*  DENS-AIR  DENSITY  AT  INTERCEPT  ALTITUDE  (KG/M3) 

C*  YIELD-THERMONUCLEAR  YIELD  (KT) 

C*  NALPHA-NEUTRON  PK  FUNCTION  ALPHA 

C*  NBETA-NEUTRON  PK  FUNCTION  BETA 

C*  CEP-INTERCEPTOR  CIRCULAR  ERROR  PROBABLE 

C*  F-NEUTRON  FLUENCE  (NEUTR0NS/CM2 ) 

C* 

C**  INPUT  DATA  ** 

CEP-600 . 

PI-3.1415927 
PERKT -3  * 16E23 

C*#  4000  FT  DENSITY  INDICATIVE  OF  MX  BASING  AREA  #* 

DENS-. 65312 
YIELD-5. 

DATA  C/- . 6775E 1 1 . 5269E-2 » - . 54364E-5 1 - . 21468E-3 1 - . 38214E1 1 
t . 1 0875E2 t - . 1 3975E 1 / 

NALPHA-34 » 539 
NBETA-2 . 242 

C**  CALCULATE  SLANT  RANGE  USING  10  CELL  MODEL  ** 

DO  200  I-1'10 

CALL  TENCELL < SR »CEPi I) 

C*»  CALCULATE  MASS  INTEGAL  ** 

IF<SR.LT . . 1 )  THEN 
SR-.001 
ENDIF 

C»*  CHANGE  SR  TO  KILOMETERS 
SR-SR*. 0003048 
MI-(DENS)*SR 
MI-MH100 . 

C**  COMPUTE  HOMOGENEOUS  AIR  4PIR2  FLUENCE  ** 

PIR2F»EXP<C(1)+C(2>*MI+C(3)*MI**2.+C(4)*MI**1.5+C(5)*MI**.5+ 
*C(6)*MI**(l./3.)+C(7)*(L0G(MI) >) 

C**  COMPUTE  FLUENCE  ** 

F-< (PIR2F)/<4.#PI*(SR**2, ) )#(PERKT)*( YIELD) *( IE-10) ) 

C**  Z-NORMAL  TABLE  ENTRY  VALUE  FOR  NEUTRON  FLUENCE  ** 

Z- ( LOG  <  F ) -NALPHA ) /NBETA 
CALL  PROB(ZfCPK) 

200  CONTINUE 

RN-DRAND<6) 

IF(RN.LE.CPK)  THEN 
CALL  PATRB(0»4) 

CPK-0 ♦ 

RETURN 

ENDIF 

CPK-O. 

C 

C**  RV  KILL  MECHANISM  ** 

100  DTIM-DTIM+. 08333 
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IF ( GATRB ( 3 ) » EQ • 3 . )  THEN 
L-L+l 
ENDIF 
C 

Ct*  RV  OVERPRESSURE  PK  ** 

CttGLOSSARYt* 

Ct  OPALPHA-OVERPRESSURE  PK  FUNCTION  ALPHA 
C*  OPBETA-OVERPRESSURE  PK  FUNCTION  BETA 

Ct  SOP-SCALED  OVERPRESSURE 

C*  OP-OVERPRESSURE 

C 

C**  INPUT  DATA  tt 
CEP-1215.2 
Y-IOOO. 

OPALPHA-6 . 8755 
OPBETA-. 12435 

Ctt  CALCULATE  SLANT  RANGE  USING  TEN  CELL  MODEL  ** 

DO  300  1-1 .10 

- CALC  TENCELC CSR. CEP  » I ) 

C**  SCALE  SR  FOR  PRESSURE  AND  YIELD  t* 

SR-SR*<  <12.69/14.7>tt<l./3. >>/< < 1 .5*Y>**< 1 ,/3. > ) 
Ctt  CHECK  FOR  CELL  1  ** 

IF  (SR.LT..1)  THEN 

S0P-1E5 

ELSE 

Ctt  CHANGE  SR  TO  KILOMETERS  ** 

SR-SR*. 0003048 

SOP-EXP  < . 1 9*  <  LOG  <  SR ) ) **2- 1 ♦ 5*L0G ( SR ) - . 1 ) 

ENDIF 

C«*  SCALE  OVERPRESSURE  FOR  ALTITUDE  ** 

OP-SOP* ( 12.69/14.7) 

C«*  Z-NORMAL  TABLE  ENTRY  VALUE  FOR  RV  OVERPRESSURE  ** 
Z- <  LOG  <  OP ) -OPALPHA ) /OPBETA 
CALL  PROB<Z»CPK) 

300  CONTINUE 
ROPK-CPK 
CPK-O. 

C 

Ct*  RV  CRATERING  PK  ** 

C**GLOSSARYt* 

C*  RA-APPARENT  CRATER  RADIUS 
C*  RE-RADIUS  OF  EJECTA 

Ct  RCPK-RV  CRATER  PROB  OF  KILL 

RA— 484.54 
RE-1041.76 
RE-REt ( 2 . /3 « ) 

IF ( GATRB ( 3 ) . EQ • 2 )  THEN 
RCPK-. 1704 

ELSE1F<GATRB<3) »EQ*3)  THEN 

RCPK-. 1971 

ENDIF 
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PK-ROPK+  < 1 -ROPK ) * ( RCPK > 

C**  COMPUTE  CORRECT  PK  FOR  OU  SHELTER  «* 

IF ( GATRB ( 3 ) . EQ . 2 )  THEN 

DUPK-PK 

RN-DRAND<&> 

C**  is  OU  DESTROYED?  ** 

IF(RN.LE.DUPK)  THEN 
DU-1 
RETURN 
ENDZF 
RETURN 
ENDIF 

C**  COMPUTE  CORRECT  PK  FOR  MX  SHELTER  ** 

IF ( GATRB  <  3 ) . EQ ♦ 3 )  THEN 

IF(L»EQ*2)  THEN 

MXPK-1- < 1-MXPK ) **2 

ELSEIF(L.EQ.l)  THEN 

MXPK-PK 

ENDIF 

RETURN 

ENDIF 

C**  IS  MX  DESTROYED?  ** 

5  IF(NTC<54) . EQ.NTC< 19) )  THEN 
JI-0 
L-0 

RN-DRAND(6) 

IF(RN.LE.MXPK)  THEN 

PRINT*? '  MX  DESTROYED' 

KK-KK+1 

IF ( NRUN .  GE . NRUNS )  THEN 
MXPS- 1 . - <  FLOAT  <  KK  > /FLOAT  <  NRUNS  > ) 

PRINT*?  '»>MX  SURVIVABILITY-  ' ?MXPS 
ENDIF 
ELSE 

IF < NRUN. GE. NRUNS)  THEN 
MXPS-1 . - <  FLOAT  <  KK ) /FLOAT  <  NRUNS ) ) 

PRINT*? '»>MX  SURVIVABILITY-  '?MXPS 
ENDIF 
ENDIF 
ENDIF 
RETURN 
END 

SUBROUTINE  UI 

COMMON/USER/MISS  ?  X  ?  Y  ?  ZZ  ?  DU  ?  DUPK  ?  MXPK  ?  JI ?  L  ?  PK 
INTEGER  DU? JI ?L 

REAL  MI8S?X?Y?ZZ? DUPK ? MXPK ? PK ? CPK ? Z 
MISS-O. 

X-Y-ZZ-O. 

DU-0 

DUPK-O. 

MXPK-0 ♦ 
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JI-L-0 

PK-CPK-Z*0. 

RETURN 

END 

SUBROUTINE  TENCELL  <SR»CEP»I> 

C 

COMMON/USER/MI SS rXrY  rZZ? DU t DUPK  »  HXPK  » JI »L»PK 
C 

CttGLOSSARYt* 

C*  THETA < I) -REFERENCE  ANGLE  TO  RADII  TO  CELL  CENTERS  MEASURED 

C*  FROM  MAJOR  AXIS  COUNTER  CLOCKWISE 

C*  RHO< I) -CELL  RADII  (RHO/CEP) 

C 

CttDECLARE  VARIABLES 
INTEGER  I 

REAL  THETA(10)fRH0(10)'CEP'SR 
Ct*  DATA  INPUT  ** 

DATA  THETA/O* 45  » 135 » 225 »315»0r72' 144 t 216 » 288/ 

DATA  RHO/O t . 7 109  > . 71 09 1 . 7109 » . 7 1 09 1 1 . 509  » 1 . 509 1 1 . 509  > 1 . 509 » 

*1.509/ 

Ct*  CALCULATE  INTERCEPTOR  SLANT  RANGE  ** 

SR-RHO(I)*CEP 

RETURN 

END 

C 

SUBROUTINE  PROB(ZrCPK) 

COMMON/USER/MISS  »X» Y»ZZ»DU» DUPK  »  MXPK  r  J I  ■»  L » PK 
C 

CttGLOSSARYt* 

C*  PK-PROBABILITY  OF  KILL 

C*  CPK-CUM  PROBABILITY  OF  KILL 

Ct*  DECLARE  VARIABLES  ** 

REAL  PKtCPK'Z 
IFtZ.LT .0. )  THEN 
Z»ABS(Z> 

Ct*  CURVE  FIT  FOR  STANDARD  NORMAL  CURVE  ** 

PK- . 5*  < 1 .  +  . 19685*Z+ . 1 151 94*Z*Z+ . 000344*Z*Z*Z+ . 019527*Z**4 . ) **  < -4 . > 
ELSE 

PK-1- .5*  < 1 .  +  . 196854*Z+ . 1 15194*Z*Z+ . 000344*Z**3 .  + . 019527*Z**4 . > 

***<-4.) 

ENDIF 

PK-PK/10. 

CPK-CPK+PK 

RETURN 

END 
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SUBROUTINE  US(ISNfDTIH) 

COMHON/USER/M I SS »  X  t Y  r  ZZ  t DU 1 1 DU2  » DUPK * MXPK » JI » L  *  PK » J J 
COMMON/QVAR/NDE t NFTBU < 1 00  >  # NREL (100), NRELP ( 1 00 ) > NREL2 < 100 ) t 
INRUN r NRUNS »NTC( 100) * P ARAN < 1 00 » 4 ) *  T  BEG  r TNOW 
INTEGER  JI'L'I' JJ»DU1 » DU2»KK 

REAL  MISS  »  ZZ  »  X  r  Y  »  NO  »  DUPK  t MXPK  »  PK  »  RA  »  RE » RCPK » XX  r MXPS 
REAL  PIR2F  f HI t PI » PERKT » SR » DENS » F  *  C  <  7 ) » Z » CPK » RATT  <  5 ) 

REAL  SOP»OP r CEP t YIELD r OPALPHA » OPBETA » NALPHA t NBETA 
DATA  CPK»KK/0.»0/ 

GO  TO  <lr2»3t4f5»6)fISN 

1  DTIM-1. 

RETURN 

C**  ENSURES  CORRECT  RV  TARGETING.  *  SHELTERS  >  1/2  •  RVS. 

2  DTIH-0.0 
TR-REHST<22) 

IF ( (NTC( 32) . EQ. 21 ) .AND. <NTC( 33>  «EQ. 1 > .AND. (NTC(35) *EQ . 1 ) ) THEN 
CALL  ST  AGO ( 22  »  34 1 TR  f  0 » RATT ) 

ELSEIF ( ( NTC ( 32 ) .EQ.21) .AND. (NTC(33) .EQ . 1 > . AND. <NTC<34> .EQ. 1 > >THEN 
CALL  ST  AGO (22f35»TR»0r  RATT ) 

ELSE IF < (NTC<32). EQ.21) .AND. (NTC(34) .EQ. 1 ) .AND. <NTC<35> .EQ. 1 > )THEN 
CALL  STAGO <22»33»TR»0» RATT ) 

ELSE IF < (NTC<32) . EQ.20) .AND. (NTC(33> .EQ.2) .AND. (NTC(35) »EQ.l) >THEN 
CALL  STAGO < 22 p 34 »TR»0» RATT) 

ELSEIF<<NTC< 32). EQ.20). AND. <NTC<35> .EQ.2) .AND. <NTC<34) .EQ. 1 ) )THEN 
CALL  ST  AGO  <  22 » 33  f  TR » 0  »  RATT ) 

ELSEIF  < <NTC<32) .EQ.20) .AND. <NTC<34) .EQ.2) .AND. <NTC<35) .EQ. 1 ) )THEN 
CALL  ST  AGO  <  22  »  33  *  TR  »  0  »  RATT  ) 

ELSEIF< <NTC<32) .EQ.20) .AND. <NTC<34) .EQ.2) .AND. <NTC<33) .EQ.l) )THEN 
CALL  STAGO < 22 » 35 »TR» Ot RATT) 

ELSEIF < (NTC (32) .EQ.20) .AND. <NTC<35) .EQ.2) .AND. <NTC<33) .EQ. 1 ) )THEN 
CALL  ST  AGO  <22r34»TR»0f  RATT ) 

ELSEIF ( (NTC (32) .EQ.20) .AND. <NTC<33) .EQ.2) . AND. <NTC<34) .EQ.l ) )THEN 
CALL  ST  AGO  <  22  r  35  >  TR » 0  »  RATT ) 

ELSEIF ( (NTC(32). EQ.l?). AND. <NTC< 33) .EQ.2) .AND. <NTC(34> .EQ.2) )THEN 
CALL  STAGO  <  22  »  35  »  TR » 0  f  RATT ) 

ELSEIF  < <NTC<32) .EQ.l?) .AND. <NTC<35) .EQ.2) .AND. <NTC<34) .EQ.2) )THEN 
CALL  ST  AGO  <  22  f  33  »  TR 1 0  »  RATT ) 

ELSEIF  < <NTC<32) .EQ. 19) .AND. <NTC<33) .EQ.2) .AND. <NTC<35) .EQ.2) )THEN 
CALL  ST  AGO  <22*34*TR»0»  RATT ) 

ENDIF 

RETURN 

3  DTIM-1. 

RETURN 

C**  SEND  RV  TO  LOADS  INTERCEPTORS  ** 

C**  IS  RV  TARGETED  AT  HPST  ** 

4  IF  < GATRB <3)*EQ.l.)  THEN 

DTIH-O. 

RETURN 

ENDIF 


DTIM-0.0 

IF< (DU2.GT . .5) .AND. <GATRB<3) . EQ.4) )  GO  TO  110 
1F( JX .GE.2)  THEN 

IF < ( GATRB ( 3 ) . EQ • 2 ) .OR* (GATRB<3> .EQ.4) )  THEN 
00  TO  110 
ENDXF 
ENDXF 

C«<  CHECK  STATUS  OF  DU1  ** 

JI-JI+1 

IF < (0U1 »GT« .5) .AND. (DU2.0T » .5) >  00  TO  110 
IF  < (DUl.GT. .5) .AND. (DU2.LT. .5) )  THEN 
CALL  NODMOD ( 53 » 55 ) 

RETURN 

ENDIF 

C**  TWO  INTERCEPTORS  FIRED  ** 

IF < ( JI.GE.2) .AND. (DU2.LT. .5) .AND. < JJ.LT.3) )  THEN 
C**  IS  THIRD  RV  TARGETED  AT  A  DU?  *« 

CALL  NODHOD ( 53  >  55  > 

ENDIF 

C**  ALL  INTERCEPTORS  LAUNCHED  ** 

1F< JI.GE.4)  THEN 
GO  TO  110 
ENDIF 

C**  INTEREPTOR  NEUTRON  PK  ** 

C**GLOSSARY*« 

C*  HI-MASS  INTEGRAL  (G/CM2) 

C*  PIR2F-4*PI*R*  2  FLUENCE 

C*  PERKT-NEUTRONS  PER  KILOTON 

C*  INT ALT-ALTITUDE  OF  INTERCEPT  <KM> 

C*  SR-DI STANCE  BET WEN  BURST  AND  TARGET  (KM) 

C*  DENS-AIR  DENSITY  AT  INTERCEPT  ALTITUDE  (KG/M3) 

C*  YIELD-THERMONUCLEAR  YIELD  (KT) 

C*  NALPHA-NEUTRON  PK  FUNCTION  ALPHA 

C*  NBET A-NEUTRON  PK  FUNCTION  BETA 

C*  CEP-INTERCEPTOR  CIRCULAR  ERROR  PROBABLE 

C*  F-NEUTRON  FLUENCE  (NEUTR0NS/CM2) 

C* 

C**  INPUT  DATA  ** 

CEP-600. 

PI-3. 1415927 
PERKT-3* 16E23 

C**  4000  FT  DENSITY  INDICATIVE  OF  MX  BASING  AREA  ** 

DEN3-. 65312 
YIELD-5. 

DATA  C/- . 6775E1 1 . 5269E-2  » - . 54364E-5  » - . 21468E-3  *  - . 38214E1 1 
I . 10875E2 t - ♦ 1 3975E 1 / 

NALPHA-34 »  539 
NBET A— 2.242 

C**  CALCULATE  SLANT  RANGE  USING  10  CELL  MODEL  ** 

DO  200  X-l»10 
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CALL  TENCELL <  SR  >  CEP » I ) 

C«*  CALCULATE  HASS  INTEGAL  ** 

IF (SR.LT * . 1 )  THEN 
SR-.001 
ENDIF 

C**  CHANGE  SR  TO  KILOMETERS 
SR-SR*. 0003048 
HX-(DENS>*SR 
MI-MI*100. 

C**  COMPUTE  HOMOGENEOUS  AIR  4PIR2  FLUENCE  ** 

PIR2F-EXP<C<l>+C<2>*MI+C<3>#MI**2.+C<4>*MI**1.5+Ct5)*MI**.5+ 

*C<6)*MI**<l./3.)+C<7)*<L0G<MI>>> 

C**  COMPUTE  FLUENCE  ** 

F«<<PIR2F)/<4.*PI*(SR**2. > >*<PERKT)*<YIELD>*< IE-10) ) 

C**  Z-NORMAL  TABLE  ENTRY  VALUE  FOR  NEUTRON  FLUENCE  ** 

Z- <  LOG  <  F ) -NALPHA  > /NBETA 
CALL  PROB(ZtCPK) 

200  CONTINUE 

RN-DRAND<6> 

IF(RN.LE.CPK)  THEN 
CALL  PATRB<0» 4) 

CPK-O. 

RETURN 

ENDIF 

CPK-O. 

C 

C**  RV  KILL  MECHANISM  ** 

110  IF ( GATRB <  3  )  • EQ ♦  3  )  THEN 

L-L+l 
ENDIF 
C 

C**  RV  OVERPRESSURE  PK  ** 

CttGLOSSARYtt 

C*  OPALPHA-OVERPRESSURE  PK  FUNCTION  ALPHA 

C*  OPBETA-OVERPRESSURE  PK  FUNCTION  BETA 

C*  SOP-SCALED  OVERPRESSURE 

C*  OP-OVERPRESSURE 

C 

C**  INPUT  DATA  ** 

CEP-1215,2 

Y-1000. 

OP ALPHA-6. 8755 
OPBETA-. 12435 

C**  CALCULATE  SLANT  RANGE  USING  TEN  CELL  MODEL  ** 

DO  300  1*1 >10 

CALL  TENCELL (SR>CEP> I) 

C**  SCALE  SR  FOR  PRESSURE  AND  YIELD  ** 

SR-SR* < < 12 .69/14.7) **< l./3.))/((l ♦ 5*Y) **< 1 ,/3 *  > ) 

C«*  CHECK  FOR  CELL  1  ** 

IF  <  SR . LT . . 1 )  THEN 
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SQP-1E3 

ELSE 

Ct*  CHANGE  SR  TO  KILOMETERS  ** 

SR-SR*. 0003048 

SOP-EXP ( . 1 9* ( LOG ( SR  >  >  **2- 1 . 5*L0G  <  SR ) - . 1 ) 

END  IF 

Ct*  SCALE  OVERPRESSURE  FOR  ALTITUDE  ** 

OP-SOP*( 12.69/14. 7) 

Ct*  Z-NORMAL  TABLE  ENTRY  VALUE  FOR  RV  OVERPRESSURE  ** 
Z- <  LOG  <  OP  > -OPALPHA ) /OPBETA 
CALL  PROB<Z»CPK> 

300  CONTINUE 
ROPK-CPK 
CPK-O. 

C 

Ct*  RV  CRATERING  PK  ** 

CttGLOSSARYt* 

Ct  RA-APPARENT  CRATER  RADIUS 
Ct  RE-RADIUS  OF  EJECTA 

C*  RCPK-RV  CRATER  PROB  OF  KILL 

RA-484.54 
RE-1041.76 
RE-RE* < 2. /3.) 

IF  < <GATRB<3) .EQ.2) .OR. <  GATRB  <  3  > .EQ.4) )THEN 
RCPK-. 1704 

EL8EIF<GATRB<3) .EQ.3)  THEN 

RCPK-.1971 

ENDIF 

PK-ROPK+  < 1 -ROPK ) *  <  RCPK ) 

Ct*  COMPUTE  CORRECT  PK  FOR  DU  SHELTER  ** 

IF ( GATRB (3). EQ.2)  THEN 

DUPK-PK 

RN-DRAND <  6  > 

Ct*  IS  DU1  DESTROYED?  ** 

IF(RN.LE.DUPK)  THEN 
DU1-1 

IF  < ( DU1 • GT » . 5 ) * AND • ( DU2 » LT »  »  5  > >  THEN 
CALL  NODMOD <  53 » 55 ) 

ENDIF 

RETURN  '  '  * 

ENDIF 

RETURN 

ENDIF 

IF < GATRB (3). EQ.4)  THEN 
DUPK-PK 
RN-DRAND<6) 

Ct*  ID  DU2  DESTROYED?  t* 

IF(RN.LE.DUPK)  THEN 
DU2-1 

IF  <  <  DU2 • GT • . 5 ) • AND .(DU1.LT..5))  THEN 
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CALL  N0DM0D<55.53) 

EHDIF 

RETURN 

ENDIF 

RETURN 

ENDIF 

C**  COMPUTE  CORRECT  PK  FOR  MX  SHELTER  ** 

IF<GATRB<3) »EQ*3)  THEN 

IF (L»EQ»2)  THEN 

MXPK-1-<1-MXPK>**2 

ELSEIF (L.EQ. 1 )  THEN 

MXPK-PK 

ENDIF 

RETURN 

ENDIF 

C**  IS  MX  DESTROYED?  ** 

5  XX-NTC<54>  +  NTC<56) 

DTIM-. 167 

IF  <NTC< 19) • EG • XX) THEN 
JI-0 
JJ-0 
L-0 

RN-DRAND ( 6 ) 

IF  <  RN » LE ♦ MXPK )  THEN 

PRINT* »'  MX  DESTROYED7 
KK-KK+1 

IF(NRUN.GE.NRUNS)  THEN 
MXPS-1 . - ( FLOAT  <  KK  > /FLOAT ( NRUNS ) > 

PRINT*. '»>MX  SURVIVABILITY-  '.MXPS 
ENDIF 
ELSE 

IF (NRUN.GE. NRUNS)  THEN 
MXPS-1 . - ( FLOAT ( KK ) /FLOAT  <  NRUNS  > ) 

PRINT*.  '»>MX  SURVIVABILITY-  '.MXPS 
ENDIF 
ENDIF 
ENDIF 
RETURN 

C**  IS  RV  TARGETED  AT  MPS?  ** 

6  IF ( GATRB  <  3 ) . EG . 1 ) THEN 

DTIM-O* 

RETURN 
ENDIF 
DTIM-O. 0 

IF < (DU1 «GT» .5) .AND. (GATRB(3> .EG .2) )  GO  TO  100 
IF(DU2.GT. .5)  GO  TO  100 
C**  HAVE  TWO  INTERCEPTORS  BEEN  FIRED?  ** 

IF  <  JJ.GE . 2>THEN 

C**  IS  THIRD  RV  TARGETED  AT  A  DU?  ** 

IF  < <GATRB<3> .EQ.2) . OR . <  GATRB (3>.EQ.4>) THEN 
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60  TO  100 
ENDIF 
ENDIF 
JJ-JJ+1 

Ct*  ALL  INTERCEPTORS  LAUNCHED  ** 

IF<  < JJ.GE.3) .AND. (DU1 .LE» .5> )THEN 
CALL  N0DM0D<55r53) 

ENDIF 

IF ( JJ.GE.4)  GO  TO  100 
Ct*  NEUTRON  KILL  FUNCTION  ** 
rat* 

CEP-400. 

PI-3.1415927 
PERKT-3.16E23 
DENS- ♦ 65312 
YIELD-5. 

DATA  C/- . 6775E1 » . 5269E-2  » - . 54364E-5  f - . 21468E-3 » - . 38214E1 » 
t ♦ 1087SE2  f- . 13975E1 / 

NALPHA-34 . 539 
NBETA-2.242 

Ct*  CALCULATE  SLANT  RANGE  USING  10  CELL  MODEL  ** 

DO  500  I-1'10 

CALL  TENCELL < SR »CEP» I > 

Ct*  CALCULATE  NASS  INTEGRAL  «** 

IF<SR.LT • . 1 )THEN 
SR- .001 
ENDIF 

Ct*  CHG  SR  TO  KM 

SR-SR*. 0003048 

MI-DENS*SR 

HI-MI*100. 

Ct*  COMPUTE  HOMOGENEOUS  AIR  4PIR2  FLUENCE  ** 

PIR2F»EXP<C<1)+C<2)*MI+C<3)*MI**2.+C<4)*MI**1.5+C<5>*MI**.5+ 
tC(6)*MI**(l./3. >+C(7)*<L0G<MI ) ) ) 

Ct*  COMPUTE  FLUENCE  ** 

F»<(PIR2F)/(4.*PI*(SR**2.))*(PERKT)*YIELD*1E-10) 

Ct*  Z-NORMAL  TABLE  ENTRY  VALUE  FOR  NEUTRON  FLUENCE  ** 

Z«  <  LOG  <  F ) -NALPHA ) /NBETA 
CALL  PROB(ZfCPK) 

500  CONTINUE 

RN-DRAND<6) 

IF ( RN . LE » CPK  >  THEN 
CALL  PATRB<0» 4) 

CPK-0 . 

RETURN 

ENDIF 

CPK-O. 

Ct*  RV  KILL  MECHANISM  tt 
100  IF(GATRB<3) .EG. 3)  THEN 
L-L+l 
ENDIF 


c**  RV  OVERPRESSURE  PK  ** 

CEP-1215.2 

Y-IOOO* 

OP ALPHA-6 . 8755 
OPBETA-. 12435 

C**  CALCULATE  SLANT  RNG  USING  10  CELL  MODEL  ** 

DO  600  1-1 >10 

CALL  TENCELL <  SR » CEP  » I  > 

C**  SCALE  SR  FOR  PRESSURE  AND  YIELD  ** 

SR-SR*<  < 12.69/14. 7  >**< l./3»>>/(<l . 5*Y>**< 1 ./3» ) > 
C**  CK  FOR  CELL  ONE  ** 

IF(SR.LT..l) THEN 
S0P-1E5 
ELSE 

C**  CHG  SR  TO  KILOMETERS  ** 

SR-SR*. 0003048 

SOP-EXP ( . 19* (LOG (SR ) )**2 . -1 «5*L0G<  SR>- » 1 ) 

ENDIF 

C**  SCALE  OVERPRESSURE  FOR  ALTITUDE  ** 

OP-SOP* <12. 69/14. 7) 

- ZSTnJGT0PT-0PALPHA770PBETA 

CALL  PROB(Z'CPK) 

600  CONTINUE 
ROPK-CPK 
CPK-0 . 

C 

C**  RV  CRATERING  PK  ** 

RA— 484.54 
RE-1041.76 
RE-RE* < 2. /3.) 

IF < (GATRB<3> .EQ.2) .OR. <GATRB(3> .EQ.4) )  THEN 
RCPK-. 1704 

ELSEIF  <  GATRB  <  3  > »EQ.3)THEN 
RCPK-. 1971 
ENDIF 

PK-ROPK+ < 1 -ROPK ) *RCPK 
C**  COMPUTE  CORRECT  PK  FOR  DU  SHELTER  ** 

IF  <  GATRB  <  3 ) »  EQ • 2 ) THEN 

DUPK-PK 

RN-DRAND<6) 

C**  IS  DU1  DESTROYED?  ** 

IF(RN.LE.DUPK)THEN 

DU1-1 

RETURN 

ENDIF 

RETURN 

ENDIF 

IF (GATRB (3) .EQ.4)  THEN 
DUPK-PK 
RN«DRAND<6> 
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C*t  IS  DU2  DESTROYED?  ** 

IF  <  RN . LE • DUPK ) THEN 
DU2-1 

IF < (DU2.GT .  «5> .AND. (DU1 »LT. .5) )  THEN 
CALL  NODHOD  <  55  f  53  > 

ENDIF 

RETURN 

ENDIF 

RETURN 

ENDIF 

C**  COMPUTE  CORRECT  PK  FOR  MX  SHELTER  ** 

IF  <GATRB<3) . EQ.3>THEN 
IF(L.EQ.2)THEN 

MXPK=l-< 1-MXPK ) **2 
ELSEIF<L.EQ» 1 )THEN 
MXPK-PK 
ENDIF 
RETURN 
ENDIF 
END 

SUBROUTINE  UI 

COMMON/USER/MISS  f  X f Y  f  ZZ » DU 1 1 DU2 f DUPK  f  MXPK  f  JI f  L  f  PK  f  J J 
INTEGER  DUfJIfLfJJ 

REAL  MI SS  f X  f  Y  f Z2  f  DUPK  f  MXPK  f PK  f  CPK  f  Z 
MISS-O. 

X-Y-ZZ-O. 

DU1-DU2-0 

DUPK-O. 

MXPK-O. 

JI-JJ-L-0 
PK— O. 

CPK=0 . 

Z-0. 

RETURN 

END 

SUBROUTINE  TENCELL  (SRfCEPfI) 

C 

COMMON/USER/MISS  f  X  f  Y  f  ZZ  f  DU1 f DU2 1 DUPK  f  MXPK  f  JI f  L  f  PK  f  J J 
C 

C«*GLOSSARY«* 

C*  THET A  < I > -REFERENCE  ANGLE  TO  RADII  TO  CELL  CENTERS  MEASURED 

Ct  FROM  MAJOR  AXIS  COUNTER  CLOCKWISE 
C*  RHQ< D-CELL  RADII  <RHO/CEP> 

C 

CMDECLARE  VARIABLES 
INTEGER  I 

REAL  THETA<10)fRH0(10)fCEPfSR 
C**  DATA  INPUT  ** 

DATA  THETA/O f 43 f 1 33 f 225 f315f0f72f!44f216f 288/ 

DATA  RHO/Of .710? f ,710? f .7109 f .710? f 1 .50?f 1 ,509f 1 ,509f 1 .50?f 
*1.30?/ 
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CS*  CALCULATE  INTERCEPTOR  SLANT  RANGE  ** 

SR-RHO(I)*CEP 

RETURN 

END 

C 

SUBROUTINE  PROB(ZrCPK) 

COMMON/USER/M I SS »X»Y»ZZ»DU1* DU2 » DUPK » MXPK » J 1 1 L  t PK , J J 
C 

CBtGLOSSARYt* 

C*  PK=PROBABILITY  OF  KILL 

C*  CPK*CUM  PROBABILITY  OF  KILL 

C«*  DECLARE  VARIABLES  ** 

REAL  PK  * CPK  * Z 
IF<Z.LT.O.)  THEN 
Z*ABS(Z> 

C**  CURVE  FIT  FOR  STANDARD  NORMAL  CURVE  ** 

PK* . 5*  < 1 .  + . 19685*Z+ . 1 151 94*Z*Z+ . 000344*Z*Z*Z+ . 01 9527*Z**4 . ) ** ( -4 . ) 
ELSE 

PK*1- .5*< 1 .+ . 196854BZ+ . 1 15194*Z*Z+ . 000344*Z**3 . + . 019527*Z**4 ♦ ) 
t**<-4.) 

ENDIF 

PK-PK/10. 

CPK»CPK+PK 

RETURN 

END 
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Figure  15-1.  Q-GERT  Network  (Layered  Defense) 


Figure  15-2.  Q-GERT  Network  (Layered  Defense) 


Figure  15-3.  Q-GERT  Network  (Layered  Defense) 


Figure  16-1.  Q-GERT  Network  (One  TDU  Defense) 


Figure  16- 


Figure  16-3.  Q-GERT  Network  (One  TDU  Defense) 


Figure  17-1.  Q-GERT  Network  (Two  TDU  Defense) 


Figure  17-2.  Q-GERT  Network  (Two  TDU  Defense) 


Figure  17-3.  Q-GERT  Network  (Two  TDU  Defense) 
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Model  Modifications 

Layered  and  One  TDU  Defense  Models.  If  the  number  of 
shelters  is  greater  than  or  equal  to  the  number  of  attacking 
RVs,  then  the  number  of  servers  for  activities  10,  11,  and  12 
must  be  changed  to  the  number  of  shelters  minus  two,  one,  and 
one,  respectively.  When  the  number  of  attacking  RVs  is  one 
more  than  the  number  of  shelters,  the  number  of  servers  should 
be  modified  as  stated  above  and  user  subroutine  two  (US2) 
should  be  changed  to  the  following: 

2  DTIM=0 . 

TR=REMST(22) 

IF  (GATRB  (4)  ,EQ.  23)  THEN 

IF ( (NTC (32) . EQ . 2 1 ) . AND . (NTC (33) . EQ. 1) .AND. (NTC (34) .EQ. 1) )THEN 
NREL2  (32)=*  1 
NREL2 (33) =1 
NREL2 (34)»1 
END  IF 
END  IF 
RETURN 

Two  TDU  Defense  Models.  If  the  number  of  shelters  is 
greater  than  or  equal  to  the  number  of  attacking  RVs,  then 
the  number  of  servers  for  activities  10,  11,  98,  and  12  must 
be  changed  to  the  number  of  shelters  minus  three,  one,  one,  and 
one,  respectively.  When  the  number  of  attacking  RVs  is  one  or 
two  more  than  the  number  of  shelters,  the  number  of  servers 
should  be  modified  as  stated  above  and  user  subroutine  two 
(US2)  should  be  changed  to  the  following: 
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2  DTIM-0.0 

TR-REMST (22) 

IF (GATRB (4) . EQ . 14)THEN 

IF ( (NTC (32) .EQ. 11) .AND. (NTC(33) .EQ.l) .AND. (NTC  (34) .EQ.l) .AND 
§  (NTC (35) . EQ. 1) )THEN 
NREL2 (32)*1 
NREL2 (33)*1 
NREL2 (34)*1 
NREL2 (35)«1 
END  IF 

ELSEIF (GATRB(4) .EQ. 15)THEN 

IF ((NTC (32) .EQ. 12) .AND. (NTC (33) .EQ.l) .AND. (NTC (34) . EQ. 1) . AND 
§  (NTC (3  5) . EQ. 1) )THEN 

NREL2(32)-1 
NREL2 (33)*1 
NREL2 (34)*1 
NREL2 (3S)*1 
END  IF 
END  IF 
RETURN 


Q-GERT  Network  Documentation 

Layered  Defense  Model. 

(1)  RV  Generation  (Node  1)  -  Generates  the  required 
number  of  RVs  using  attribute  1  as  a  counter.  These  RVs 
arrive  normally  distributed  in  time  with  a  mean  of  five  seconds 
and  a  standard  deviation  of  two  seconds. 

(2)  Attribute  Assignment  (Node  1)  -  Attributes  2  and  4 
are  set  equal  to  the  number  of  shelters  minus  two  and  the 
number  of  shelters,  respectively.  Both  attributes  are  used 
for  probabilistic  branching  in  the  model. 

(3)  Take  First  Branching  (Node  3)  -  If  the  number  of 
RVs  is  twice  the  number  of  shelters,  the  RVs  branch  to  node 
41  for  random  assignment  to  a  shelter.  If  the  number  of  RVs 
is  less  than  twice  the  number  of  shelters,  the  RVs  branch 
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to  node  30  to  perform  an  equivalent  procedure.  The  number 
of  RVs  will  never  be  greater  than  twice  the  number  of  shelters 
due  to  the  fratracide  limit  previously  mentioned. 

(4)  RV  Target  Assignment  (RVs  Twice  the  Number  of 
Shelters  -  At  node  41,  the  RVs  are  probabilistically  targeted 
(using  attributes  2  and  4)  to  either  an  empty  MX  shelter  (node 
42),  the  terminal  defense  unit  (node  43;  located  in  one  of  the 
shelters),  or  the  actual  MX  (node  44),  and  are  assigned  an 
attribute  3  of  1,  2,  or  3,  respectively.  The  number  of  servers 
and  queue  capacities  are  restricted  and  balkers  are  returned 

to  node  41  to  insure  that  two  RVs  are  targeted  to  each  shelter. 
From  here,  all  RVs  are  sent  to  the  interceptor  queue  (node  19) 
with  a  flight  time  delay  of  ten  minutes. 

(5)  RV  Target  Assignment  (RVs  Less  Than  Twice  the  Number 
of  Shelters)  -  At  node  30,  the  RVs  are  sent  to  node  31  with  a 
time  delay  of  one  minute.  At  node  31,  the  RVs  are  probabilisti¬ 
cally  targeted  (using  attributes  2  and  4)  to  either  an  empty  MX 
shelter  (node  32),  a  terminal  defense  unit  (node  33),  or  the 
actual  MX  (node  34),  and  assigned  an  attribute  3  of  1,  2,  or  3, 
respectively.  When  an  RV  is  assigned  to  a  target,  it  encounters 
user  subroutine  two  (US2)  which  insures  that  a  minimum  of  one 
RV  is  targeted  at  each  shelter.  The  number  of  servers  and 
queue  capacities  are  restricted  and  balkers  are  sent  back  to 
node  31  which  insures  that  all  protective  shelters  are  targeted. 
After  this  procedure,  all  RVs  are  routed  to  the  interceptor 
queue  (node  19)  with  a  flight  time  delay  equal  to  ten  minutes. 
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(6)  Set  RV  Status  -  At  queue  node  19,  attribute  4  is 
reset  equal  to  one  to  indicate  that  all  RVs  are  active. 

(7)  Select  RVs  to  be  Intercepted  by  Exoatmospheric  Layer  - 
At  node  20,  if  an  RVs  attribute  1  is  less  than  or  equal  to 

the  number  of  interceptors  available,  the  RV  will  be  sent  to 
queue  node  21  to  await  interaction  with  an  exoatmospheric 
interceptor.  Although  the  RVs  are  intercepted  first -come- 
first -serve,  the  exoatmospheric  defense  does  not  know  the  RV's 
intended  target  since  the  targets  were  randomly  assigned  in 
steps  4  and  5  above.  If  an  RV  does  not  encounter  an  exoatmos¬ 
pheric  interceptor,  the  RV  will  penetrate  this  layer  of  defense 
and  await  interception  by  the  endoatmospheric  layer  which 
consists  of  one  terminal  defense  unit  (TDU). 

(8)  RV/Exoatmospheric  Interceptor  Engagement  -  Each  RV 
in  queue  node  21  is  serviced  by  an  exoatmospheric  interceptor 
in  accordance  with  user  subroutine  three  (US3),  which  deter¬ 
mines  if  the  RV  is  destroyed  using  a  CEP  distribution  coupled 
with  a  "cookie-cutter"  weapon  radius  as  mentioned  in  the 
structural  model. 

(9)  Determine  RV  Status  -  After  being  serviced  by  the 
exoatmospheric  layer,  each  RV  is  examined  by  user  subroutine 
one  (US1)  which  routes  those  RVs  not  destroyed  to  node  51  to 
join  those  RVs  which  penetrated  the  exoatmospheric  defense  and 
are  awaiting  interaction  with  the  TDU.  From  node  51,  all  RVs 
are  sent  to  the  TDU  queue  (node  53)  with  a  flight  time  delay 
of  three  minutes. 
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(10)  RV/Endoatmospheric  Interceptor  Engagement  -  Each 
RV  in  queue  node  53  is  serviced  by  an  endoatmospheric  inter¬ 
ceptor  in  accordance  with  user  subroutine  four  (US4)  which 
employs  the  strategy  stated  in  the  structural  model.  US4 
determines  if  an  RV  has  destroyed  the  TDU  and  the  RV’s  proba¬ 
bility  of  killing  (PK)  the  MX.  If  the  TDU  is  destroyed,  the 
MX  cannot  be  defended  by  this  layer  of  defense. 

(11)  MX  Status  -  After  all  of  the  RVs  have  been  serviced 
by  the  defense,  user  subroutine  five  (US5)  compares  the  proba¬ 
bility  of  killing  the  MX  to  a  random  number  to  determine  if 
the  MX  is  destroyed.  US5  computes  the  number  of  times  the  MX 
is  destroyed  in  N  simulation  runs  and  prints  out  MX  surviva¬ 
bility. 

Endoatmospheric  Defense  Model  (One  TDU). 

(1)  through  (6)  -  Same  as  the  layered  defense  model. 

(7)  Flight  Time  -  All  Rvs  in  queue  node  19  are  sent  to 
queue  node  53  with  a  flight  time  delay  of  three  minutes. 

(8)  The  remaining  steps  are  the  same  as  steps  (10) 
through  (11)  of  the  layered  defense  model. 

Endoatmospheric  Defense  Model  (Two  TDUs). 

(1)  Same  as  the  layered  defense  model. 

(2)  Attribute  Assignment  (Node  1)  -  Attributes  2  and  4 
are  set  equal  to  the  number  of  shelters  minus  three  and  the 
number  of  shelters,  respectively.  Both  attributes  are  used 
for  probabilistic  branching  in  the  model. 
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(3)  Same  as  the  layered  defense  model. 

(4)  RV  Target  Assignment  (RVs  Twice  the  Number  of 
Shelters  -  At  node  41,  the  RVs  are  probabilistically  targeted 
(using  attributes  2  and  4)  to  either  an  empty  MX  shelter 
(node  42),  the  primary  TDU  (node  43),  the  backup  TDU  (node  45), 
or  the  actual  MX  (node  44),  and  are  assigned  an  attribute  3  of 
1,  2,  4,  or  3,  respectively.  The  number  of  servers  and  queue 
capacities  are  restricted  and  balkers  are  returned  to  node  41 
to  insure  that  two  RVs  are  targeted  to  each  shelter.  From 
here,  all  RVs  are  sent  to  the  interceptor  queue  (node  19)  with 
a  flight  time  delay  of  ten  minutes. 

(5)  RV  Target  Assignment  (RVs  Less  Than  Twice  the  Number 
of  Shelters  -  At  node  30,  the  RVs  are  sent  to  node  31  with  a 
time  delay  of  one  minute.  At  node  31,  the  RVs  are  probabilisti¬ 
cally  targeted  (using  attributes  2  and  4)  to  either  an  empty 

MX  shelter  (node  32),  the  primary  TDU  (node  33),  the  backup 
TDU  (node  35),  or  the  actual  MX  (node  34),  and  are  assigned  an 
attribute  3  of  1,  2,  4,  or  3,  respectively.  When  an  RV  is 
assigned  to  a  target,  it  encounters  US2  which  insures  that  a 
minimum  of  one  RV  is  targeted  at  each  shelter.  The  number  of 
servers  and  queue  capacities  are  restricted  and  balkers  are 
sent  back  to  node  31  which  insures  that  all  protective  shelters 
are  targeted.  After  this  procedure,  all  RVs  are  routed  to 
the  interception  queue  (node  19)  with  a  flight  time  delay  equal 
to  ten  minutes. 
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(6)  Same  as  the  layered  defense  model. 

(7)  Flight  Time  -  All  RVs  in  queue  node  19  are  sent  to 
node  53  with  a  flight  time  delay  of  three  minutes. 

(8)  RY/ Endoatmospheric  Interceptor  Engagement  -  Each  RV 
in  node  S3  is  routed  to  node  54  with  no  time  delay.  Each  RV 
in  node  53  is  serviced  by  the  primary  TDU  according  to  US4 
which  utilizes  the  strategy  outlined  in  the  simulation  section 
of  this  study.  If  the  primary  TDU  is  destroyed  or  has  launched 
two  of  its  three  interceptors,  node  53  is  replaced  by  node  55, 
and  the  backup  TDU  assumes  the  defense  role  if  it  has  not  been 
destroyed.  If  both  TDUs  are  destroyed,  the  MX  cannot  be 
defended  by  this  defense  structure.  The  backup  TDU  services 
each  RV  according  to  user  subroutine  six  (US6).  If  the  backup 
TDU  is  destroyed  or  has  launched  all  three  of  its  interceptors, 
node  55  is  replaced  by  node  53  and  the  primary  TDU  reassumes 
the  defense  role  if  it  has  not  been  destroyed. 

(9)  Same  as  step  (11)  of  the  layered  defense  model. 
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Appendix  H 


List  of  Acronyms 

Antiballistic  Missile 
Antisatellite 
Ballistic  Missile  Defense 
Circular  Error  Probable 
Intercontinental  Ballistic  Missile 
Infrared 

Intensity  Sure -Safe 
Intensity  Sure-Kill 
Kilotons 

Low  Altitude  Defense 

Multiple  Protective  Shelters 

Missile  X 

MX  Survivability 

Neutrons  per  square  centimeter 

Probability  of  Kill 

Pounds  per  square  inch 

Reentry  Vehicle 

Terminal  Defense  Unit 

Weapon  Radius 
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